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I. INTRODUCTION 

In this MS thesis, we have followed Eric Gourgoulhon's lecture notes titled 3+1 Formalism and Bases 
of Numerical Relativity (arXiv:gr-qc/0703035vl lUl and E. Poisson's book : A Relativist's Toolkit, The 
Mathematics of Black-Hole Mechanics, Cambridge University Press, Cambridge (2004) fH. 

Historically, the 3+1 approach has been put forward G. Darmaois (1927)[5|, A. Lichnerowicz (1930-40) 
JU, Q and Y. Choquet-Bruhat (1952)18]. During 1958, 3+1 formalism started to be used to construct the 
Hamiltonian form of general relativity by P. A. M. Dirac 111, lITOl and later by R. Arnowitt, S. Deser and C. 
W. Misner (1962) flF]. The 3+1 formalism became popular in the numerical relativity community during 
1970 ID. 

The 3+1 formalism is used to rewrite the Einstein equation as an initial value problem and construct 
the Hamiltonian form of the general relativity. This method is based on the concept of the hypersurface, Z,, 
which is independent of whether the given spacetime is a solution of the Einstein equation or not. In this 
formalism, we consider that there is an embedding mapping O which maps the points of a hypersurface into 
the corresponding points of the four-dimensional manifold M such that M is covered by the continuous set 
of hypersurfaces (J^t)teR- Furthermore, the well-known Gauss-Codazzi relations and the 3+1 decomposition 
of the spacetime Ricci scalar curvature are the fundamental equations of the 3+1 decompositions of the 
spacetime (Al, g). And they play a crucial role in the 3+1 decompositions of the Einstein equation. The 
Gauss-Codazzi relations are defined on a single hypersurface. On the other hand, the 3+1 decomposition 
of the spacetime Ricci scalar is obtained from the flow of the hypersurfaces. Moreover, the foliation is 
valid for any spacetime with a Lorentzian metric so we have to restrict our selves to the globally-hyperbolic 
spacetimes. And the foliation kinematics of the globally-hyperbolic spacetimes allow us to construct the 
Ricci equation whose contraction with respect to the induced 3-metric gives the last fundamental equation 
of the 3+1 formalism (i.e the 3+1 expression of the spacetime scalar curvature) HI. 

The 3+1 decomposition of the Einstein equation is obtained by using the Gauss-Codazzi relations, the 
3+1 decomposition of spacetime Ricci scalar and the 3+1 decomposition of stress-energy tensor. Basi- 
cally, the four-dimensional Einstein equation decomposes into three main equations which are known as: 
the dynamical Einstein equation, Hamiltonian constraint and Momentum constraint. The dynamical Ein- 
stein equation is obtained from the full projection of the Einstein equation onto the hypersurface and has 6 
independent components, the Hamiltonian constraint is obtained from the full projection of Einstein equa- 
tion along the normal vector and has 1 independent component and the Momentum constraint is obtained 
from the mixed projection of the Einstein equation and it has 3 independent components. Therefore, as we 
expect, the total number of independent components are 10 which is exactly the number of independent 
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components of the Einstein equation in four-dimensional spacetime |[T1. 

Tiie 3+1 dimensional Einstein system is modified to the Cauchy problem (or initial- value problem) 
by rewriting it as a set of PDEs (Partial Differential Equations) and specifying with the help of particular 
choices of the lapse function N and shift vector (5. Choosing a scalar field N, a vector field yS and a spatial 
coordinate system (x') on an initial hypersurface allows us to define a unique coordinate system {x") within a 
neighborhood of Eq such that x° = corresponds S,. That is , N and yS are depend on the coordinate systems. 
And also, the lapse function N at each point of Eq leads us to define a unique vector m(= N n) which is 
used to construct the neighboring hypersurface l^st by Lie dragging each point of So along m. Therefore, 
the 3-1-1 dimensional Einstein system can be turned into as a PDEs system by using tensor components 
which are expanded with respect to the coordinates {x") - (t, x') adapted to the foliation. The PDEs form 
of the 3-1-1 dimensional Einstein system contains only tensor fields of Sf and their time derivatives which 
implies that they can be taken as a time evolving tensor fields on a given Ej. The PDEs form of the 3-1-1 
dimensional Einstein system is a system of second-order, non-linear PDEs for unknown (yij , Kij , N , fi' ) 
when the matter source terms {E , pi , S ij ) are given. Here y/y is the metric of the hypersurface, Kij is the 
extrinsic curvature, E is the energy density, pi is the momentum density and S ij is the stress tensor. The 
crucial point is that the PDEs form of the 3-1-1 Einstein system contains neither the time derivative of N nor 
of p. This means that they are not dynamical variables rather they are just quantities associated with the 
coordinates {x") = {t,x') (that is, Lagrange multipliers). Therefore, PDEs form of the 3-1-1 Einstein system 
can be converted into the initial value problem by choosing particular N and yS HI, H. 

Beside the orthogonal decomposition that is used for 3+1 formalism, the conformal decomposition is 
also used to define the flow of the hypersurfaces by continuously mapping an initial well-defined conformal 
background metric jij into the induced 3-metric of X,. Lichnerowicz |l6l proposed that by mean of the 
particular conformal decomposition of the extrinsic curvature, one can arrange the constraint equations 
which allows us to define initial data for the Cauchy problem. In addition to this, York has shown that 
the conformal decompositions can be used for the time evolution |[T5l . That's, he has proved that the two 
degrees of freedom of the gravitational field are carried by the conformal equivalence classes of the induced 
3-metric lITSll . The Weyl tensor is used to check whether a given spacetime, whose dimension is greater than 
3, is conformally flat or not. And it disappears for lower-dimensional manifolds. In this case, the Cotton- 
York tensor { II14L II15L II13I1 }. C'^, does the same task of the Weyl tensor in higher dimensional spacetime. 
Furthermore, the Cotton- York 11141 . ifTSl . |[T3]| tensor of weight 5/6, - y^/^C'-' , is conformally invariant. 

The Hamiltonian model approaches a physical state at a certain time and gives the evolution of the state 
as time varies. This model is being transformed into the gravitational theory as a state on a particular 
spacelike hypersurface ||9]|, ifTOll . Now, the gravitational theory is a covariant theory and locally has Lorentz 
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symmetry. In order to write the Einstein equations into the Hamiltonian form, people started to give up the 
main covariance property of the gravitational theory by choosing a family of particular coordinate systems 
such that " jc" = constant" corresponds a spacelike hypersurface. Instead of the set (yij , Kjj , N , fi' ) in the 
PDEs form of 3+1 Einstein system, Arnowitt, Deser and Misner have proposed the ADM formalism of the 
general relativity in which conjugate momentum of the induced three-metric jij, tt'^ = ^/y{Ky'-' - K'-'), is 
used. In the ADM formalism |[T6) . n'-i and yij are the dynamical variables and the Lapse function N and the 
shift vector yS are taken as Lagrange multipliers |[T6ll . In the chapter 6, we will see this in detail iQl. Pll. 

The action for the General Relativity (when the boundary term is different than zero) contains Einstein- 
Hilbert part and Matter part. The infinitesimal four-dimensional volume element is taken as the union of two 
spacelike hypersurfaces Z,, , which are at the upper and lower boundaries and a timelike hypersurface S 
between Z,, , E,^. Now, 3+1 decomposition of M and 2+1 decomposition of the timelike hypersurface with 
proper choice of vectors lead us to the conserved quantities of the ADM mass, ADM linear momentum and 
ADM angular (by using rotational Killing vectors) of a given hypersurface. However, due to the fact that 
global quantities of mass, linear momentum and angular momentum are defined only for asymptotically 
flat spacetimes, the ADM formulas are vahd just for the spacetimes which asymptotically converge to well- 
defined space times such as Minkowski spacetime HI, Q, II16I . 

Finally, R. Bartnik and R Tod introduce ||20l the conditions on the intrinsic quantities of the Zf. They 
ensure whether the is a hypersurface of a spacetime which is a solution of the four-dimensional static 



vacuum field ( with A = 0) or not. In addition to [20], we introduce the equations (404i, (405 1 (for A = 



case) and (402 1, (403 1 (for KtO case). These equations can be used to find which solutions of the gradient 
Cotton (Conformal) soliton [21] are also the solutions of the static vacuum fields equations. However, We 
have not also been able to solve the constraint equations and have not found explicit metric. Moreover, we 
think that our method can be extended to the Ricci soliton ||23]| and Topologically Massive Gravity (TMG) 
m, l2ll,l|26l. 



II. GEOMETRY OF HYPERSURFACES 

Since 3+1 decomposition of the spacetime is constructed by slicing the spacetime with a continuous 
set of the hypersurfaces, {Lt)teR, we will deal with the geometrical fundamentals of the hypersurface. The 
geometrical results that we will obtain in this chapter are fully independent of whether the given spacetime 
is a solution of the Einstein equation or not. The only constraint is that the spacetime must have Lorentzian 
metric El, S, . 
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A. Notations and Basic Geometrical Tools 

1. Spacetime and Tensor Fields 

We assume a real, smooth {i.e. C°°) four-dimensional manifold M which endowed a Lorentzian metric of 
signature (-, +, +, +) and a connection V. In general, no one can define a global vector space on manifolds. 
Therefore, it is considered that at each point of the manifolds there is a space of vectors Tp(M) (titled 
as tangent space at the point p) and corresponding space of linear forms T*(M) (titled as dual space or 
cotangent space at the point p ). Furthermore, we suppose that all the Greek letters [a,p,y,...} run in 
{0, 1,2,3} are free indices , {yu,y,p, ...} are dummy indices and all Latin letters [i,j,k, ...} run in {1,2, 3} and 
{a, b, c, ...) run in |2, 3). 

Since the Tp{M) and T*{M) are vector spaces, we consider that there is a set of basis {e") which spans 
Tp{M) and the associated dual set of basis {ca) which spans T*{M) such that e"{ep) - 5" p. Therefore, any 
tensor field T of type can be expanded with respect to these bases as 

T = T"'-"<'p,...p^ea, ® ... ® ea^ ® ® ... ® A . (1) 

Here T^^-^p are the related components of T relative to the bases {e°') and (£„). A tensor field T 
with rank (^j turns into another tensor field VT with the rank {^^^^ when the covariant derivative acts on it. 
Therefore, the expansion of VT in these bases is 

VT = r'^'-%,...^,;ye„i ® ... ® ea^ (8 e^' ® ... ® A ® . (2) 

The contraction of covariant derivative of the tensor field T with an arbitrary vector field u gives us the 
covariant derivative of T along the vector field u which does not change the rank of tensor fields T 

V„T = VT(...,...,u), 

(3) 

p + q slots 

where u^'^ nT"^-"p p^_p^ are the components of VuT with respect to {e") and {ea). 

2. Scalar Products and Metric Duality 

In general, we do not know how to relate the elements of Tp{M) ( or of T*{M) ). The concept of metric 
is introduced to do this task. Now, the scalar product of two vectors is taken place by mean of the related 
metric g of the manifold At 

^{vi,y)eTp{M)®Tp{M), 
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u.v = giuf'ei.yey) = uf'v^gie^^ey) = g^yU^'v" = m^v^ . 

Here the metric g is taken as if it has two slots for inputting vectors. Alternatively, the same job is done by 
bracket notation : V(w, v) e T*{M) » Tp{M), 

< W, V > WfjC'^, V^Cy > 

= < e^,ey > 

= WpV^S'^y 

= Wf,vf . 

As we see in equation (j2]), " V^Wo-e" ® " are the components of the 2-form Vw relative to the bases e"^ 
of T*{A\)®T*{A\). Then, the directional covariant derivative of a 1-form iv along a vector field m, Vuiv, is 
a 1-form 

VuW = Vw(u) 

= [VyWjse'^e^Yu^'ey) 

Since the directional derivative, VyW, is a 1-form, we use the bracket notation to get a scalar from it: 
V(w, u, v) e T*{M) ® T{M) ® T{M) , 

Vw(u,v) VuW, V > 

=< u^V^wpe^, v^Cy > 
= u^VpWpv^ <e^,ey> 

= U^V^V pWy . 

Any element of Tp{M) (or T*{M)) can be mapped into T*{M) (or Tp{M)) by mean of the 2-form g. 
That's, the metric g induces an isometry between Tp{M) and T*{M). Some of crucial properties of this 
isometry are 

1. The dual of any vector u (e Tp{M)) is a unique linear form of Tp{M) and denoted by u such that 
the scalar product is defined as 

V V € Tp{M) , < u, V >^ g(u, v) . (4) 
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2. The dual of any linear form w ( € 7~p(At)), w, is a unique vector w e 7~p(At) such that 

Vv € Tp{M) , g(w, v) =< w, V > . 

3. T : TpiM) ® Tp{M) — > */? (i.e any rank tensor T maps 2- vectors of tangent spaces at the point p 
into the space of scalars). An endomorphism T is induced from T such that T : T{M) T{M) and 
it satisfies 

T(U,V) - u"\/^Ty^da^6j3'^ 

^u"v''Ty^<e\ea> (5) 
= u%''T^^g{ey,ea) 

- gyau^v'^ry^ 

= u.T(v). 

As we see in the equation ([l| and because the endomorphism T is a 1-form, T"p are the components 

—* 

of T relative to {ca)- 



3. Curvature Tensor 



According to gravity, the matter curves the geometry and the geometry determines the motion of the 
matter. The rank Riemann curvature tensor measures how much the spacetime is curved. Basically, it is 
a map which sends a 1-form and 3 vectors into the real, smooth space of scalar fields C°°{M,K) 

: T*{M) ® T{Mf C°°{M, !R) , 



(w, w, u, — >< w, VuVyW - VyVuW - V[u v]W > . (6) 

The Riemann tensor is assumed to be machine which has 1 slot for 1-form and 3 slots for vectors. The 
relation (j6|l is nothing but a tensor field on At. Furthermore, '^Wsap is the components of ^R with respect to 
a proper set of basis {Ca) and {e") of Tp{M) and T*{M). Now, the crucial properties of "^R are 

1 . "^Rafjys ( - gay'^Ky i}y6 ■ ) is anti-symmctric between the 1 ^* two terms a and jS and between the 2""^ two 
terms y and 6. 
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2. ^Rafjys satisfies tlie cyclic property between tlie last three indices which is known as second Bianchi 
identity 

^RajiyS + ^RaySp + ^RaSPy - . (7) 

3. For the torsion-free spacetime, the well-known relation of the Ricci identity is 

[V„,V^]w^-^/?VyS>v^- (8) 

4. One-times contraction of^Ryasp , i.e. 5 — > yS, leads us to a symmetric, bilinear-form Ricci tensor ^R. 
The '^R is considered to be a machine that has 2 slots for vectors 

^R : TiM) ® T{M) ^"(M, K) , 

(9) 

(u,v)^4^(e^u,ep,v). 
Also, the trace of the Ricci tensor relative to the dual of g results in the spacetime Scalar curvature R. 

5. The traceless part of the spacetime Riemann tensor titled as Weyl tensor, which gives whether 
a given spacetime is conformally flat or not is obtained by subtracting all the trace part (i.e. the 
Ricci tensor) and the trace-trace part (i.e. the Ricci scalar = g'^^R^iv ■ ) of the spacetime Riemann 
curvature tensor from itself 

'^C^SaP =^R^Sap - ^{^R^agSp " ^R^ pgSa + ^RspS^ a " '^RdaS^ jj) 

1 ^ (10) 

+ ■^^R[g6aS^li - gSpS^a) ■ 

We need to emphasize that the Weyl tensor , "^C, vanishes for spacetime whose dimension is lower 
than 4. Thus, in the lower dimensional geometry the spacetime Riemann tensor can be written in 
terms of the Ricci tensor, the metric and the scalar curvature tensor 

4. Hypersurface Embedded in Spacetime 

As we see in the figure ([T]), the set of points, Vp € M, at which the scalar field is constant corresponds 
a hypersurface X of the four-dimensional manifold At which is an image of a three-dimensional manifold 
S under the homeomorphism <t. Since the three-dimensional manifold E is something like to be embedded 
into the four-dimensional manifold At, we say that the mapping O is an embedding mapping 

^■.t^M. (11) 

Furthermore, "one-to-one character of the embedding mapping O ensures that the hypersurfaces do not 
intersect" . 
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FIG. 1 . The embedding of t, into M 



The embedding mapping O induces two well-known mappings of the push-forward mapping, , and 
the pull-back mapping, O* . First, O* maps the vectors of the tangent space of the three-dimensional mani- 
fold X, Tp(%), into the corresponding vectors of Tp{M) 

<D. : Tp{t) Tp{M) , 

(12) 

V = (v^, v^ v^) (Dh-v - (0, v^ v^) . 

Here V - (v^, v^, v^) is the components of the vector v with respect to the natural basis d/dx'of 7~p(At) 
associated with the coordinates (x'). On the other hand, O* maps the linear forms of TpiM) into the 
corresponding linear forms of Tp(%) 

<:>* ■.T;iM)^T;it), 

w 0*w : Tpit) n , (13) 

Further insight, the mapping acts on the multilinear forms of Tp{M), too 

V(vi, v„) e , {^*T){yy, v„) = r(<D,vi, O.vJ , (14) 

where T is an n-form. Especially, the pull-back of the 2-form spacetime metric g takes a great attention. The 
pull-back of g is called the induced metric, y, of the hypersurface, Ef and known as the first fundamental 
form of X, 

r = a)*g. (15) 
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Moreover, the scalar product between any two vectors of the tangent space of the hypersurface is same 
either by using g or y: V(u, v) e TpCS) x 7~p(S) , 

u.v = g(u, v) = 7(u, v) . (16) 
5. Normal Vector 

We consider there is a scalar field tonM such that each "t ^constant" corresponds the hypersurface Sf 
of M and the vector Vt is normal to S. Then, the dual of the vector Vt is the gradient 1-form dt such that 
the relation between them is 

V^t^g'^^V^t^g'^'^idt^. (17) 

and V V e X, the scalar product vanishes 

<d/,v>=0. (18) 

Furthermore, the type of the normal vector Vt is determined by the type of the hypersurface E: That's, if the 
signature of induced metric y of X is (+,+,+) then E is spacelike and the corresponding normal vector Vt 
is timelike, contrary, if the signature of induced is (-,+,+) then S is timelike and the corresponding normal 
vector Vt is spacelike and, finally, if the induced metric is degenerate , i.e. has signature of (0,+,+) then 
either T,t or Vt are null. 

Although Vt is a unique normal vector to S,, it is not a unit normal vector. Therefore, in the not-null 
case, we normalize it to get a unit normal vector of n 

Vt 

n - , , (19) 

■\/±Vt.Vt 

where the positive sign (+) is used for a timelike hypersurface and the negative sign (-) is used for a spaceUke 
hypersurface. Thus, the norm of the unit normal vector is 

Vt.Vt |~1 if li is spacelike 
n.n = ±- 



Vt.Vt 



+1 if 1, is timelike 



6. Intrinsic Curvature 



For the not-null case, one can always propose a unique Levi-Civita connection D which is still torsion- 
free and metric compatible on a hypersurface S. Moreover, the intrinsic covariant derivative D is defined 
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by using the induced metric y. Now, as we see in the relation ([6]l, the Riemann curvature tensor measures 
the curvature by using the spacetime connection V. However, in order to measure the curvature of the 
hypersurface Z (i.e. the intrinsic curvature of the E), we replace the spacetime connection V with the 
intrinsic connection D. 

V V e rCI) , [D,-, Dj\v^ - R^ujv^ ■ (20) 

That's, the intrinsic curvature of a given hypersurface is nothing but the curvature which is measured (or 
felt) by anybody moving on the hypersurface. Also, as we did in the equation (|9]l, one-times contraction on 
the intrinsic curvature tensor gives us the Ricci tensor of the hypersurface S 

^ij - R^ikj ■ (21) 

Finally, the contraction of the intrinsic Ricci tensor results into the intrinsic scalar curvature (or Gaussian 
curvature) of S. 

R = y'jRij. (22) 



7. Extrinsic Curvature 

In 3+1 formalism, the global manifold is assumed to be constructed by a family of embedded hypersur- 
faces. Naturally, we expect that there must be a machine (or a tensor field) which will measure how much 
the hypersurfaces are bending within the global manifold. Fortunately, there is one which is known as the 
extrinsic curvature K. To find the explicit form of K let us first image a physical case: when a drop of ink is 
being released orthogonally onto the surface of water, it spreads over the surface. Therefore, the orthogonal 
release of the drop is something related to the spread of it over the surface. By taking this approach as a 
reference, we see that we need a vector which is related to the spread of the unit normal vector n over the 
hypersurface E. Basically, this is done by endomorphism Weingarten map (or shape operator) x of Tp(L) 

X:Tp{i:)^Tp(L), 
V — > Vvfi . 

In words, the machine x inserts the unit normal vector n into its slot and migrates the directional derivative 
ofn which is an element ofTpCL) 

n.;^(v)-n.V,n- iv,[n.n] -0. 
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Now, let us deduce the crucial property "self-adjointness with respect to the induced metric y " of x'- 

v(u,v)erp(X)®rp(i), 



Uykf(v) - u-Vyft - Vv[u.n] - n.VyU - -n.VyU , 



We assume that the torsion tensor is zero (VyU - VuV - [u, v] = O). Then, ( 

U-Ar(v) = -n.(VuV- [u,v]) 

- -n.VuV + n.[u, v] 

= -Vu(n.v) + v.VuA + n.[u, v] 

- v.Vun + n.[u, v] . 
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) becomes 



For the sake of self-adjointness of x, we need to show that the last term of ( 24 1 disappears 



V?.[u,v] ^< dt,\u,y] > 

^< <if,(VuV - VyU) > 

dt, VuV > - < dt, VyU) > 

^< V^te^',u''VyV^es > - < V ^.t e'' yVyU^eg > 
^ Vf.tu^VyV^ <e>',es> -Vf^tv^VyU^ < e'',es > 



= v^tu^vysi's - v^.tv^vys^'s 

= Wf,tu''VyV^ -V^tV^'WyU'' 

= 0, 

where we used the fact that is orthogonal to V^f. Thus, we proved that;^ is really self-adjoint 

V(u, V) € TpiY) ® r/I) , n.x{y) = Ar(u).v , 



(23) 



(24) 



(25) 



(26) 



Since the Weingarten map^ is self-adjoint, its eigenvalues are taken as the principal curvatures, Ki, of the 
hypersurface S and the corresponding eigenvectors are taken as the principal directions of the hypersurface 
such that the mean of the Ki is known as the mean curvature,'H, of X 

1 ^ 

•^ = -2^,. (27) 

i=i 

Contrary to the intrinsic curvature, /c, and 'H are depend on how the hypersurface is embedded into M so 
they are taken as extrinsic character of X. Now, we are ready to construct the explicit mapping of 2-fomi of 
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the extrinsic curvature K of S by using the Weingarten map;^f. K is assumed to be a machine which has two 
slots for the vectors of Tp(L) and whose output is a scalar 

K:7-p(i:)®7-/I)^^, 

(28) 

(u, v) — > -n.x(y) . 

This is the well-known relation of the second fundamental form (or the extrinsic curvature tensor ) of the 

hypersurface E. Moreover, the relation between the contraction of K with respect to the induced metric y 
and the mean curvature, of E is 

K - f^Kij - -S-K . (29) 

Up to now, we dealt with timelike and spacelike hypersurfaces. Now, we will restrict ourself to spacelike 
hypersurface in which the signature of the induced metric y is (+,+,+) and we will define the fundamental 
geometrical tools for it: 



8. The Orthogonal Projector 

The tangent space of A1 at a point p, 7~p(A1), can be orthogonally decomposed into the corresponding 
tangent space of the hypersurface Z at the point p, 7~p(X), and a one-dimensional vector space of n, 'yect(n) 

Tp{M) = Tp(L) ® 'Vectin) . (30) 

where 'yect(n) is a 1-dimensional vector space for n. Because in the null case 'yect(n) c T"p(E), the 



orthogonal decomposition of vector space ( 30 1 is valid only for spacelike and timelike hypersurfaces. Now, 
the orthogonal decomposition ( 30 1 of T"p(Al) allows us to define an operator y which projects the elements 



of Tp(M) intoofTp(I) 

y : TpiM) Tp{T) 

(31) 

V — > V -I- (n.v)n , 

here y is known as the orthogonal projection operator. It selects the components of the vector of S among 
of At. Therefore, the projection of the unit normal vector n is equal to zero [i.e. since n.n = -1, then, 
7(n) = n -I- (n.n)n = 0] and it acts as an identity operator for vectors of Tp(Z) [/.e.Vv e Tp(L) , y{y) = 
V -I- (n. v)n = v] . Further insight, the orthogonal projection operator y can be expanded relative a set of bases 
{Ca) of Tp{M) and the corresponding components are 

y"p-5''p + n"ni,. (32) 
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We mentioned in (12 1 and ( 13 1 the embedding <1) induces the push-forward mapping, and the pull-back 
mapping, O* in the given direction and does not imply in the reverse directions. On the other hand, as we 
illustrated in ( [jlj l that y carries the elements from Tp{M) and projects them into of TpCL). And, it induces 
another mapping yj^ between the corresponding dual spaces [from T*{1,) to T*{M)] 

w € , V € TpiM), y*M ■■ r;{Y.) t;{M) , 

V — >< ryviw,v > (33) 

= w(r(v)) . 

— 

Also, the induced mapping can map arbitrary n-form ^ of Tpil.) 

V j?l e 7-;(I)" ,y*M-^^ y*M^ ■ T-piM)" n , 

(vi, v„) y*M^{yi, v„) , (34) 
= J?l(y(vi),...,y(vJ). 

Particularly, the extension of 2-form induced metric y to At will act on the vectors of Tp{M) .Then, we 

— 

denote it with the same symbol, y = y^y. The relation between the extended induced metric y and the 
spacetime metric g is 

— 

y = g + n^n^ y j^y . (35) 

where fl is a 1-form. As we did before, yap = gap + nanp are the components of the extended induced metric 

— 

y {= yyviT) relative to a proper family of basis {e") of Tp{M). Let us take a look at the action of y on the 
particular cases: 

1. V(v, u) e S, then, the induced metric y and the spacetime metric g will do the same job on these 
vectors 

)^ 

y j^y{\i,y) = g(u,v)-i- < n, u >< n,v > 

= g(u,v) (36) 



2. On the other hand, if one of these vectors (consider u = An) is coUinear with n, then, the action of y 
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will be zero. That's, for any v e Tp{M) 



y(u,v) = 7;k7(u,v) 

= Ag{n, v) + /I < ft, n >< ft, V > 

= ^{g(n,v)-<ft,v>) (37) 

= /}{ < ft, V > - < ft, V > 

= 0. 



By observing the equation (32 1 and the components of the extended metric y (= y y ), we see that the 
orthogonal projection operator y is obtained from the extended metric by raising its 1 index. Indeed, we 
use the same symbol for the extension of the extrinsic curvature K to M, too: 

K=y*MK. (38) 

Finally, with the help of the orthogonal projection operator y any rank-(p+q) tensor T [ € T(M)'''SiT*{M)'^ 

— ** 

] can be converted into another tensor, y^T, of same type which is still an element of 7'(M)'' ® T'*(At)^. 
The transformation between their components is by 

G*MTr-''W.A, = y"\. -r"^v/V. -r%7'^'-"%>...v, - (39) 



9. Relation Between "K" and Vfi 



Up to now we have not said anything about the diffusion of the unit normal vector n within the neigh- 
borhood of a point p of the hypersurface. We only assumed n to be at points of the hypersurface. Basically, 
we will see that deviation of n leads us to the relation between the extrinsic curvature K and the covariant 
derivative along n. Therefore, we need to define the acceleration vector a of ft 

a = Vflft. (40) 

If we assume that ft is a 4-velocity of an observer (since ft is a timelike vector), then, ft is 4-acceleration of 
the observer. Furthermore, this deviation is an element of 7~p(Z), 

(,.a = (,.V.6=lv,»« = 0. (41) 

Now, we have emphasized that K can be taken as a machine which has two slots for inputting vectors. 
Therefore, let us extend it to M and insert two vectors belonging the tangent space of M at p: V(u, v) e 
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{yMK){u,y) = K(yiu),y{y)) 



= -y(u).V^ n 

y(v) 

^ -[u + (n.u)n].V[v+(A.v)fl]n 

= -u.Vvii - (n.u)n.Vvn - (n.v)(u.a) - (n.v)(n.u)(n.a) 

- -u.Vyii - (n.v)(u.a) 

= - < Vn(..., v)), u > - < a, u >< ft, V > , 



(42) 



where we used u.v = g(u, v) = g(u'^e^, v^Cy) = gf^yu'^v^ = u'^'v^. We know that the equation! 42 1 is valid for 
any pair of tangent vectors of M so we can drop the vectors to get the compact form as 



K = -Vft - a ® ft ^> Vfl = -K - a ® ft , 



(43) 



here the symbol (~) means dual of the vector. The components of the tensorial equation (43 1 with respect to 
a given basis of T'p(AI) are 



(44) 



The equation ( 44 1 is defined onto the four-dimensional manifold M. Then, let us pull-back it to 7'*{1,) 

7^ aj" j}^ vn^i = -y'' aj" pK^v - y'' aj" pa^riy , (45) 



notice that the last term on the right hand side of the equation (45 1 vanishes because the projection of «y 
onto 2 is zero so we get 



or in compact form 



K - 7 K = -y Vn. 



(46) 



(47) 



As we see in the equation (47l, the projected form ofVh onto E (i.e. —K) is symmetric though the four- 



dimensional Vn is not. Now, it is easy to show that the trace of the equation ( 44 1 with respect to is 

K = -l^y^ , (48) 

or in compact form 

K = -V.n . (49) 



The equation ( 49 1 gives the relation between the scalar extrinsic curvature and the divergence of the unit 
normal vector. 
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10. Relation between Connections of the Spacetime and of the Hypersurface 

Consider a tensor field T onto S. Botli tlie tensor field T and its covariant derivative DT are tensor fields 
of At, too. As we implies before, we are able to convert DT into another vector field of At which is denoted 
by 7*VT 

DT = r*VT , (50) 

or more precisely, 

DT-/Vr>. (51) 
And the related transformation of the components is given by 

^p-* I3\...j3q — Y HI---Y HpY j8i---/ PqY p^cr^ vi-.v, • 



The crucial properties of the transformation of ( 5 1 



are 



1- y ^TmT is a torsion-free connection on E, it satisfies all the defining properties of a connection 
[linearity, reduction to gradient for a scalar field, commutes with contraction and Leibniz rule]. 

2- 7 ^7mT metric compatible 

= y^ay^pY'y^plgtJiy + «p«y] ^^^^ 
= y''ay''l3'fy[^pgpv + (Vp?l^)?lv + «^i(Vp?ly)} 

= 0, 

where we have used the metric compatibility of g and the fact that the orthogonal projection of n 
onto Z vanishes. 

Finally, let us derive the relation between the connections: 

V(u, V) e TpiT) ® TpCL) , [o^yf - u'^D^v" 

- u^[6% + n-«^)v,v'^ (54) 
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Now, the last term of the equation (54 1 is nothing but the definition of the extrinsic curvature (28 1. By using 



this fact, the equation ( 54 1 becomes 



DuV = VuV + K(u, v)n . 



(55) 



That's, the difference between the directional covariant derivatives D„v and V„v is given by the extrinsic 
curvature K(u,v). Furthermore, the difference is along the unit normal vector h. 

B. Geometry of Foliation 



From now, as we see in the figure ([2]), we will deal with a continuous family of embedded hypersurfaces 
(^ijteK which covers the spacetime (yVt, g). However, the spacetime (At,g) that we are going to foliate (i.e. 
to slice) is not any type of spacetime rather it is a globally-hyperbolic spacetime. 




t+dt 



FIG. 2. The foliation of (M, g) 



1. Globally Hyperbolic Spacetimes and Foliation 

In order to define what the globally-hyperbolic spacetime is, we need to first define the concept of 
Cauchy Surface: Any spacelike hypersurface which is being intersected by causal curves (i.e. timelike and 
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null) if and only if one time is called A Cauchy Surface. Now, any spacetime on which there can be defined 
a continuous family of Cauchy surface is called a globally-hyperbolic spacetime [figure Q]. Also, it is 
obvious that the topology of the globally-hyperbolic spacetime is essentially E x In words, the foliation 
of the spacetime is that there is assumed to be a smooth scalar field t on the four-dimensional manifold At 
such that the union of points at which the scalar field is identical construct the related hypersurface 

Wte'R,'Lt = [peM,t{p) = t]^M = ^2^, (56) 

teR 

where we assume that the gradient of the scalar field is always different than zero (i.e. regular) which 
ensures that the hypersurfaces never intersect 

I.t n S^' = for ttt . (57) 



2. Foliation Kinematics 



1 . Lapse Function 

For convention, let us use the symbol t as the scalar field on M. Then, the vector V? is essentially 
normal to the hypersurface and not necessarily a unit normal vector. Therefore, we suppose that 
there exists a scalar field N{> and called lapse function ifTTTl ) which is used to re-normalize Vt : 

Vt 1 
n = -NWt = — , ^ - , , (58) 



the minus sign in the equation (58 1 guarantees that n is future-oriented. Because V? is a vector, its 
dual is the gradient 1-form dt. The crucial point is that the lapse function makes the dt a unit 1-form, 
too: Suppose that there is a scalar field 5 such that n = S dt so what is S ?, 

< n,n > =< Sdt,-NVt > 

^SN{- <dt,Vt>] 

1 (59) 
= S.N— 

]\f2 

= -1, 

so S = -N. Thus, 



n ^ -Ndt . 



(60) 
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m(q) 




FIG. 3. The Lie Dragging of v € T(E,) along m such that XmV e 7~(E,).Here the difFeomorphism O mappings the 
elements of 7~p(2,) into of TqCE,) 



2. Normal Evolution Vector 

Now, when we evaluate the inner product between the vector and the normal unit vector n, we 
will see that it is not equal to 1. This means that the normal unit vector h can not follow the flow 
of the scalar field though it locally does. Then, the modification of the n to the evolution of the 
hypersurfaces seems as a primary condition. Therefore, we propose a new normal vector m (known 
as normal evolution vector) 



m = A'^n 



m m - -A^-^ , 



(61) 



such that the normal evolution vector m is being adapted to the flow of the scalar field t, 

<dt,m> -< dt,Nn > 
= N < dt,n > 

<dt, -NVt > 
- -N^ < dt, Vt > 



(62) 



= 1 



1 



This modification provides the evolution of the hypersurfaces. That's, all the points of the initial 
hypersurface are being carried along the vector 6tm such that the union of carried points construct 
the neighbor hypersurface E^+^j [see figure |4|. This evolution of the hypersurface is known as the Lie 
dragging of hypersurfaces along m. As we see from the figure Q the Lie dragging along m does 
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not disturb the elements tangent to S,. That is , they are still the elements tangent to after dragging 
along m: 

vv € , x^v € . (63) 




3. Eulerian Observers 

Actually, the unit normal vector n can be taken as the 4-velocity of the Eulerian Observer. Then, the 
world line of Eulerian Observer will intersect the hypersurface Z, only one time which says that all 
the simultaneous events at each constant scalar field (i.e. local) construct the hypersurface Z,. Now, 
assume that there are two points (=events) p (€ Zj) whose coordinate time is t and p (€ whose 
time coordinate i?,t + 5t on the world-line of a Eulerian Observer. Then, the elapsed time difference 
6t between two events with respect to the observer's clock is given by 

5t - yj-giSt m, 6t m) - Vg(-nMn)<5f, (64) 



from the equation (61 1, we get 

5T = N6t. (65) 



As in the equation(40 1, the corresponding 4-acceleration of the the Eulerian Observer is an element 
of TpCEt) 

a-Vftfi. (66) 



Due to the fact that a € Tp{l,t), the 4-acceleration vector a can be rewritten in terms of the lapse 



25 



function N: 



- -n'^V^m-^nal-Nn'^VaV^t 
= ^n^ni^iV^N) - Nn^'Vai-^n^) 



^nan^{Vf,N) - N^n'^nf.VaN + N^r 
1 ., _ 1 



or in compact form 



a-DlnA^. (67) 



Gradient of the 1-form n and m 

We are going to deduce very important two relations of the gradient of h and of m in terms of the 



extrinsic curvature and the lapse function. Firstly, the substitution of the equation (67 1 into the 



equation (43 1 gives us the gradient of the 1-form h , 

Vn ^ -K - DlnN^ ft =^ V/^na ^ -Ka/j - Da InNnp . (68) 



Secondly, the equations (61 1 and (68 1 are used to find the gradient of the dual of the m, 

Vm = V(A^n) 

= A^V(n) + n®VN 
= n{k - D In ® n) + ft (g) VA^ 
= - K - -DA^ ® ft) + ft (g) VA^ , 
so we get the gradient of the dual of m as 

Vm = -K - DA^ ® ft + ft (g) VA^ . (69) 

And its vector form as 

Vm = -K - DA^ (g) n + n ® VA^ ^ V/j/m„ ^ - 'NKafi - DaNn^ + n^VpN . (70) 
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Evolution of the Induced 3-Metric 

Under the flow of the hypersurfaces, the induced 3-metric is being carried by the Lie derivative along 



m. With the help of the equation the equation (70 1, the Lie derivative of y along m is deduced from 



so we get 



Now, 



+ 7ai,[-NK''i3 - D^Nnfi + ni'VpN] 

+ Jtipn^'VaN - Nyaf.K''^ - ya^D^Nnf^ + yaf,n''V/jN 

- Nn^'{y^na)n|i + Niian^iy^np) - NKpa - DpNiia - NK^p 
- DaNnp 

= Naati/s + NapUa - INKap - DpNua - DaNnp 

- NDa In Nnp + NDp In Nua - INKap - DpNua - DaNnp 

- DaNnp + DpNua - INKap - DpNria - DaNnp , 

Lr^Jap = -2NK,p =^ £r«y = -INK . 

■Lmjap = LNnJap = Nfi^'^ ^.y^p + y^pVaiNn'') + ya^.'^piNn^) 
= Nn'^Vf.yap + yf.pNV^ni' + y^pn^'V^N + ya^NVpn'' 

+ yai.n^'VpN 
= Nin^'^^yap + y^ip^an^ + yaf^^pn^) 

= NLaJap ■ 



(71) 



(72) 



Thus, the equations ( [71] ) and ( [72| ) lead us to rewrite the extrinsic curvature in terms of Lie derivative 
of the induced 3-metric along the unit normal vector: 



-CnTap - —Lmjap ^> Kap - --£n7ap ■ 



(73) 



Evolution of the Orthogonal Projector 

In order to find the evolution of the tensor field T of we need first to show what happens to 

the orthogonal projection operator under the flow of the hypersurfaces. Again, the evolution of the 
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orthogonal projection operator is done by the Lie derivative along m 

Because the projection of the normal vector and its dual onto the hypersurface is zero, we have 

£^y"p = N{n>'V^n")np + Nin^'V ^,np)rf + NK"p + {D"N)np 
- n^DpN - NK"p - {D"N)np 
= Na"n/} + Na/jii" - n"{DpN) - {D"N)nii 
= N^D^Nnn + N^DpNn" - n"{DfjN) - {D"N)nii 
= 0, 



or in compact form 



(74) 



The equation (74 1 implies that if initially a tensor field T is an element of hypersurface, it will remain 
to be an element of the hypersurface throughout the evolution. Moreover, let us show this in detail: 
if a tensor field T is a tensor field belonging to the hypersurface, then, the orthogonal projection 
operator acts as an identity operator. 



jT - T. 



(75) 



Now, we are seeking after the Lie dragging whether there is any projected part of the carried tensor 
field ofT along the normal unit vector h or not. If there is any, then, the Lie dragging will disturb the 
tensor T which is tangent to E to be still tangent to X during the flow of the hypersurfaces. Therefore, 
let us show on a simple sample of T of type Then, the transformation in terms of the components 
is 



Let us carry the transformation along m. 



Hence, the property of (74 1 allows to write (77 1 as 



(76) 



(77) 



(78) 
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or in compact form 

r*XmT = £™T. (79) 

Thus, the Lie derivative along m of any tensor field T tangent to is a tensor field tangent to E;. 

3. The Foliation of Hypersurfaces 

1 . Foliation of Spacelike Hypersurface X, 

The coordinates adapted to the foliation - x^iy") is assumed to be the parametrized curves where 
the parametrization is the induced coordinates (y") of the hypersurface H. The projection operator 
can be written as 

dx" 
" dy^ ■ 

Since the extended line element of 



(80) 



dsl = g,pdx"dxP = g,^[—dy"){—dy') - habdy"dy'' , (81) 
so the induced 3 -metric is 

Kb = gapeaell , (82) 

where hat is the metric component of the spacelike hypersurface Z,. Therefore, contravariant form 
of the spacetime metric can be decomposed as 

g'^/? = + , (83) 

where na is timelike normal unit vector to E,. Further inside, the spacelike hypersurface can be 
decomposed into its spacelike unit normal vector plus its boundary 2-surface Sf- suppose that the 
induced coordinates y" on is parametrized curves where the parametrization is the coordinates 
j^C^) of the 5/. Then, the corresponding projection operator can be taken as 

< = ^- (84) 

The extended line element of S t 

dsj^ = habdx^dJ" = hab(—d9^)(—de^) = aABdO^dO^ , (85) 
so the induced 2-metric is 

o-AB = habe^eB ■ (86) 



(88) 
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And the decomposition of the spacelike induced 3-metric is 

Therefore, the contravariant form of the spacetime metric can be written in terms of the timelike and 
spacelike unit normal vectors as 

= -renl3 + r"r^ + CT^\eyi)(4el) 

Because e" is tangent to and n" is normal to Ef, we have 

ra = =^ = raey" = ^ r„ ± , (89) 

where n" is the unit timelike vector which is normal to E, and is the unit spacelike vector normal 
to the boundary of E, (that is, 5 , ). 

2. Foliation of the Timelike Hypersurface S 

Let z' be the coordinates on the timelike hypersurface S H and the corresponding projection operator 
toEf 

The induced 3-metric of S is obtained from the projection of the spacetime metric g: 

Jij-gafielJ'.. (91) 
And the contravariant form of the spacetime metric can be decomposed as 

^aH ^ ^a^ ^ yj^a^^ _ (92) 

Because z' is arbitrary, let us choose them as z' = {t,e^) 

Here 9^ {A = 1,2) are the adapted coordinates of the boundary of the spacelike hypersurface Ef. 
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Then, the extended Une element of S is 

ds^ - gapdx"dx!^ 

= gai}(Nrfdt + elde^)(Nnf^dt + (/^de^) 

= gai3[N^n"nPdt^ + Nn"/j^dt + d9^ + Nn/^e^^dt + dO^ + e^ePj^de^dO^] (94) 

- -N^-dt^ + CTABdG^de'^ 

- jijdz'dz-' ■ 

Thus, we get the 2+1 decomposition of the 3-metric of the timelike hypersurface S as 

Yijdz'dz^ = -N^dt^ + o-ABdG^dO^ . (95) 



C. Coordinate Adapted to the Foliation 

Here, it is assumed that there is a coordinate system {x' = x^,x^, x^) ) on each Xf such that it smoothly 
varies during the flow of the hypersurfaces. Then, we take [x" - (t, x^)] as a proper coordinate system on 
At. 




FIG. 5. The Shift Vector j8 

Naturally, these adapted coordinate systems induce another set of coordinate systems for the T'p(Al). 
Basically, the corresponding partial derivatives are often chosen as the bases of Tp{M) : 

where (i= 1,2,3). Moreover, because of the shift vector that we are going to define, dt does not have to be a 
timelike vector. 
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1. The Shift Vector p 

In addition to the TpiM), the coordinates adapted to the foliation also induces a set of basis gradient 
1-form A-jf for T*{M) such that it obeys 

<dx",dp>=5%. (97) 

Because of < dt, dt >- 1 , 5, drags the hypersurfaces as m does, too. However, in general, they do not have 
to be collinear. They are collinear only if (x') of Z, are orthogonal to each other. Otherwise, as in the figure 
(|5]), there is assumed to be a deviation vector (shift vector jS) lilLl between them 

d, = m + ]3. (98) 

And it is easy to show that the shift vector p is an element of the hypersurface E,. As it is said before dt 



does not have to be timelike. This property is determined by the square of the /J: From the equation (98 1, 
the norm of 5, is 

^,2t = -N^ + yS.yS , (99) 

so the if yS.yS < A'^^ then dt is a timelike vector, or if p.jS > then dt is a spacehke vector and finally, if 
yS .yS = A'^^ then dt is a null vector. 



2. 3+1 Form of the Metric 

The 3+1 decomposition of the spacetime metric and the 2+1 decomposition of the induced 3-metric 
of the timelike hypersurface S play a crucial role during the construction of the Hamiltonian form of the 
general relativity. Therefore, we need to find the 3+1 form of the spacetime metric g, too: After defining the 
suitable coordinate systems adapted to M, let us expand the induced 3-metric y relative to these coordinates 
(xO of If 

y = yijdJ®dxK (100) 

Because the shift vector /3 is tangent to X,, we can raise its indices with the help of the components of the 
induced 3-metric jij in this coordinate system, 

pi^jijpL (101) 
Also, the expansion of the spacetime metric g in the corresponding coordinates is 

g = gaji dx" ® . 
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Basically, g can be assumed as a machine which has two-slots for the corresponding vectors. Therefore, the 
components of it are obtained by 



gap ^ g(da, dp). 



(102) 



Therefore, with the help of the equations (98 1 and ( 102 1, the components gap in this coordinate systems are 

gm = g(dt,t) = 'dt^t = -A^' + ^i/S' , (103) 



goi ^ g(5r, di)^ dfdi ^ im + /3).di 
-< Pjdx-', di >= Pj < dx^^di > 
-Pi, 



(104) 



/ ^ 

.?00 gOj 






. giO gij , 




Pi ytj , 



gij = gidi, dj) = di.d j = Jij. 
Thus, the 3+1 decomposition of the gap in matrix representation is 

gap 

Furthermore, let us evaluate the explicit 3+1 form of spacetime metric: 

gfivdx^ ^dx" - + I50^dt ®dt+ Pjdt » dx^ + fSidt ® dx' + jijdx' ® dx^ 

= ( - A^^ + ®dt + jijP'dt (g> dx^ + jijP'dt (g> dx' 

+ jijdx' ® dx^ , 

we get 

gfxvdx'' ® dx" = -N'^dt ®dt + yij\dx' + /?' Jf] g) \dx^ + /3^<if] . 



(105) 



(106) 



The ability of writing g in terms of its parts! 106 1 is something like the beginning of the 3+1 formalism and 
the Hamiltonian form of general relativity. Now, let us find the matrix form of the components of the dual 
spacetime metric g"^. Suppose that the dual metric has the form of 









( k \ 

a V 




gij 




yj yk 

\ ) 
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The matrix multiplication (In general, AjjBjk - Cik) between the covariant and the contravariant forms of 
the spacetime metric is given by 

Pi 7U 

The multiplication of the 1*' row of ga/j with the V column of g"'^ gives 

a/ii + Jijv^ = ^ api = -v; , (107) 



^ 


( k ^ 




1 




\ ) 




.0 



(-A^2 + pfi^)a + hjv^ = 1 ^ i-N^ + Pj/3j)a - a/3jpj = 1 



(108) 



From the equations ilOlh and (llOVll, we get a - and = The multiplication between 2"^ row of 



ga/j with the 2'"^ column of g"'^ will give us the components b^*^: 



Pi0l 



(109) 



Let us multiply ( 109 1 with y 



7''7ijb^' - 7"6ik - 7 



J — y 



Ar2 ' 



I -^i, k ^ j . 

With the help of the previous change of indices, we get 



b'J = 



(110) 



Thus, the decomposition of the components of the dual metric in matrix form is 



^ ,00 Oy ^ 



y y ^2 ; 



Notice that gij = yjj but that in general g'j + y'j. Alternatively, the dual spacetime metric g"^ in matrix 
form can be obtained from the Cramer's rule. Now, let us denote the determinant of y and g as 



g ^ detigafi) , 



(111) 



7 ^ detiyij) . 



(112) 
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Observe that because of the lapse function N and the shift vector /?, the related determinants of yij and ga/3 
are coordinate dependent. Because of this, they are not tensors rather they are tensor densities. Now, let us 
evaluate component by using the Cramer's rule, 

r = ;^ = ^ = ^. (113) 



The ([113| must be equal to the component of the dual spacetime metric that we have obtained from the 



matrix multiplication, 

8 8 

Thus, the relation between the determinant of the induced 3 -metric and of the spacetime metric is 

V=^-A^Vr- (115) 



III. THE GAUSS-CODAZZI RELATIONS AND THE 3+1 DECOMPOSITION OF SPACETIME RICCI 

SCALAR 

In this chapter we will deduce the fundamental relations that are in the center of the 3+1 formalism of 
general relativity yj. From now, unless it is emphasized, the hypersurface that we will work on must be 
taken as the spacelike hypersurface Z, (i.e. whose unit normal vector n is timelike ): 



A. Gauss and Codazzi Relations 



1. Gauss Relations 



In order to find the first fundamental equation of the 3+1 formalism (i.e. Gauss relation), we start with 



equation ( 20 1 (or the intrinsic Ricci identity) and use the related transformation ( 52 1 between connections 



of D and V: The Ricci identity on Z is 



(116) 



here v e Tp{M). Since the second term in the equation ( 1 16 1 is obtained by interchanging the indices of a 



and /3 of the first term, it is better to work just on the I term of the equation ( 1 16 1: With the help of the 
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equation ( 32 1, the 1" term of the ( 1 16 1 is obtained as follow: 

- /ar>^p(V;.«")«v/iV^v'^ + r^r>^«"(Vp«v)/iV^v'^ 
+ y'a//^r^pr"v(V^«0«iV^v'' + /„r>^pr"v«^(V^«i)V^v'* 

Because the projection of the unit normal vector (and its dual) on the hypersurface is zero, then, the previous 
equation becomes 

+ r^//^r^pr"v/iV^v^v'^ 
= -^a/?r^i«"v^v'' + K^a7''/sv''^<rnA + /.r"/^r^^v^v^v^ 



where we used the equation (46 1. Let us interchange the dummy index of ^ with A, 
D^Dpyy = -K^py^irf^y - Ky^Kp^ + y^rY" /sy^ a^ <rv' . 



(117) 



In order to find the 2'"' of the equation (|116|), it is enough to just interchange the indices a and /? of the 



equation (1171: 



(118) 
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If we interchange the indices yu and cr of ( 1 18 1 and then subtract it from the equation (1171, we will get 



The last term is nothing but the projection of spacetime curvature tensor. Then, we have 



(119) 



We can rewrite v'' as = jP^v^ . Then, the equation ( 1 19 1 becomes 



A^p, p^cr, cr^v, ^ ^ S A ^ 6 

by changing of the dummy indices, the previous equation turns into 



(120) 



since the vector v'' is arbitrary, then, we can drop it to get the full projection of the spacetime Riemann 
curvature tensor onto the hypersurface as 



(121) 



The equation ( 121 1 is known as the Gauss relation. Let us continue by contracting the indices a and y of 



Gauss relation ( 121 1: 



rpy'/sy^s ^R'af.v = Rsp + KKps - k,sK% . 



(122) 



The left hand side of the equation ( 122 1 can be rewritten in terms of the sum of the full projection of the 
spacetime Ricci tensor onto E and the mixed projection of the spacetime Riemann curvature tensor ^R: 

= Ypy^S 'Ray + yfyff0% %a^y 

= fpy% 'Rav + n'^nPfpy^s 'R%y^ 

cr <r^ p,p V 



fpy''s'Ray + y,sY'pn''n''RP 



vpcr ■ 



(123) 



The substitution of the equation ( 123 1 into the equation ( 122 1 gives 

y^py^s 'R<Ty + y^^sVpn^n' ^R^yp^ = Rsp + KK^s - K,sK"p . 
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By applying the given change of the dummy indices in the previous equation, we get the well-known relation 
of contracted Gauss relation as: 



„o"„v 4 



(124) 



Now, let us multiply the equation ( 124) by y"l^ (i.e. taking trace with the dual induced 3-metric y"^) 



it becomes 



(125) 



We need to first modify the V and 2"^ terms of the left hand side of the equation ([125 1 



y"" %y - [g"" + n^n'] ^Rpv = ^R + n^n' %y , 



(126) 



and 



vpcr 



Because the Riemann tensor is antisymmetric in its first two and second two indices, the last term in the 
previous equation is zero. Therefore, it turns into 



Y'.n'ff^R^p^ = n''n''^Ry^. 



(127) 



By substituting the related results of ( 126 1 and ( 127 1 into the main equation of (125 1, we get the another 



well-known relation of scalar Gauss relation as 



+ l^Ry^n^'n'^ =R + K^- KjjK'^ . 



(128) 



2. Codazzi Relation 



In order to find the second fundamental relation of 3-1-1 formalism (i.e. Codazzi relation), we start with 
the Ricci identity of the four-dimensional spacetime for the unit normal vector n. And, then, we will project 
it 3-times onto the hypersurface : Now, the Ricci identity in four dimension is 



(V„V^-V^V„)«^ = 4/?r^„^„-. 



(129) 
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Let us start by projecting the equation ( 129 1 three-times onto the hypersurface 



(130) 



Again, as we did in the part of Gauss relation, let us just work on the l"' term of the equation ( 130\ because 
the 2'"^ term is obtained simply by interchanging the indices a and jS of the 1 term: Therefore, 

= J^'aYpyW-^.KPy - {^^,aP)n, - aP{^,,n,)} 
= -rafpy^p^.KP, - yPafpy\{^paP)n, 
-taffifpaP{^^,n,). 

Since the projection of the dual of the unit vector is zero, the 2"'^ term on the last part of the previous 
equation vanishes. Then, it becomes 

= -D^K^p - Yaf0\aP{l,,n,) , 



(131) 



where we used the general transformation relation ( 52 1 between connections. Because of the equation (46 1, 



the equation (1311 turns into 



y\rYfiy\^v^vnP = -DaK^p + Kafi 



(132) 



As we said before the interchange of the indices a and yS of the projected equation ( 132 1 of the 1" term of 
( 129 1 gives the projected version of the 2"'^ term of the equation {129), 



(133) 



Finally, the subtraction of ( 133 1 from ( 132 1 provides us the famous Codazzi-Mainardi relation: 

f'afpy''p ^RPa^^vrf = DpK\ - D^K^p . (134) 



Now, Let us contract the Codazzi-Mainardi relation ( 134 1 on the indices a and y 

yPpfp ^RP^^yn'^ = DpK - D^K^p , 
[6Pp + nPnp]fp ^RP^^yn'' - DpK - DaK"p , 
fp^P ^R^a^^vrf + nPnpn'^fp ^RP^^y = DpK - D^K^p . 
Because the multiplication between symmetric and antisymmetric tensors are zero, the second term on the 
left hand side of the previous equation vanishes and we get 



y'pSPp 'RP^pyn" = DpK - D„K"p p 'R^yn'' = DpK - D,K"p 



(135) 



cr ^ V, y 



a ^ p., 
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and with the given change of the dummy indices in the equation ( 135 1, we reach 



(136) 



which is known as the contracted Codazzi relation. 

As we deduced above, the Gauss and Codazzi (or Codazzi-Mainardi) relations which are obtained from 
the various numbers of the projection of the spacetime Riemann curvature tensor onto or normal to a 
single hypersurface Hf That's, the full-projection of "^R onto the hypersurface E, gives the Gauss relation . 
And the three-times projections of^R onto the with the one-times projection of it along the unit normal 
vector n gives the Codazzi-Mainardi relation. 



Contrary to the Gauss-Codazzi relations, the conti^action of the Ricci equation that we will deduce in this 
section results in the third fundamental relation of the 3+1 decomposition of the spacetime Ricci scalar and 
is essentially based on the flow of the hypersurfaces . Moreover, we will see that the two-times projection 
of onto with the two-times projection of it along n will lead us to the Ricci equation. As we did in the 
previous section, the starting point is the Ricci identity in four dimension of n but, here, we will project it 
two-times onto Ef and one-times along n: The related four-dimensional Ricci identity for n is 



B. Last Fundamental Relations of the 3+1 Decomposition 



1. The Ricci Equation 



(137) 



Let us project the equation ( 137 1 two-times onto and one-times along h, 



+ V^[Z>^ InNny]} 



(138) 



= yaMn'^y^/si-^vK"^ - {^yir In N)n^ - IT lnA^(V,?i^) 



Vo-Kf^y (V^iy In N)ny IT In N{Va-ny)} 



= r«prV{-«"v,/:^^ + v,(z>' inA^) + «"v^^^ 



+ mnNn"'Vo-ny}, 



where we have used the equation (68 1 of the gradient of n. Because of the relation 



n 



Kf"^ = 0^ n^^yK^"^ = -K'^^^yn'^ , 
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the equation ( 138 1 becomes 

+ {D'' In N){Dy InN)} 
- ra^//j{^^.r[-^"v - D'' In Nn,] + V^ilT InN) 
+ rfV^Kf'y + {If lnA^)(DvlnA^)) . 



(139) 



Since the projection of dual of the normal vector is zero, the equation ( 139 1 reduces to 

-r{T^\nN){Dy \nN)\ 

- -Jai^YfiK^'aK^y + yai.ffi^viP' \^N) 

+ JapiY fin"^ <tK% + jai^fniiy' In N){Dy InN) 
= -K.^K'^p + Dp{Da InN) + /,//j«"V^/:^, 



(140) 



+ (D„lnA^)(D;3lnA^) 
= -Ka^K'^f, - l^{DpN){D,N) + ^DpDaN 

+ y'afprf^^rK^v + ^{D„N){DpN) 
= -Ka^K'^p + ^DpD^N + yf^ay'/sn'^yaK^v ■ 
The aim is to find such a projection of which is totally composed of the intrinsic quantities of S,. 
Therefore, we need to get rid off the last term of the equation ( |140 i: 

= Nn^'V^Ka/3 + Kf./il-NKi'a - ly'Nn^ + n^'VaN] 
+ Ka^l-NK''^ - iTNnp + nf^VpN] 

- Nn^'V^Ka/} - NKf./^K'', - K^fiD^Nria + K^n^^V^N 

- NK^^K^'p - K^^D^^Nnp + K„f,n>'VpN 



(141) 



= Nn'^V^.Ka/s - INKaf^K^fi - K^pE^Nua - Ka^D^Nnp , 
where we have used the fact that the projection of the extrinsic curvature along the unit normal vector is 



zero. Due to the property of (79 1, the full projection of the equation ( 141 1 onto Zf is 

'y*-CmKa/j - -CmKa/i 

- NraV/in^V^K^y - INyay'/iK^^rK'^y 

- r^f/iK^vD^Nn^ - Y'^faK^^D'^Nny 

- NraVfin^V^K^y - INrafpK^aK'^y , 



(142) 
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hence, we obtain 



(143) 



By substituting the equation ( |143[ ) into the main equation ( |140[ ), we obtain the crucial Ricci equation as 

(144) 



where we benefited from the fact that the intrinsic connection D is torsion-free (That's. DaDpf = D^Daf ; 



here, f is a scalar field). Furthermore, we can replace the first term of the equation i\T44j that is a projection 

pvcr 



of the spacetime Riemann tensor ^R'^pva- with the full projection of spacetime Ricci tensor ^R^v onto by 



using the contracted Gauss relation (124V. Now, the contracted Gauss relation that we have found is 



then 



ja^^fpn^'n'' ^R^p^ = R,p - y^y^ '^pv + KK,p - K^^Kf^p , 



(145) 



by interchanging the indices v and p of the previous equation ( 145 1 and then substituting it into the Ricci 



Equation ( 144 1, we get 



/a//jXp = -^^m^./^ - ]^D,DpN + Rap + KKap - IK^^^K^ p . 



(146) 



For convention, let us do the operation of p ^ v in the first term of the equation ( |146| ). Then, it becomes 

(147) 

This is the equation of the full projection of the Ricci tensor onto the hypersurface X,. And the compact 



y^ayp '^Rfiv - -J^£.iaKap - —DaDpN + Rap + KKap - IKa^K^ p . 



form of the equation {147) is given by 



^* 4jj ^ --£^Yi - -DDA^ + R + /:K - 2K.K 



1 



1 



(148) 



2. 3+1 Expression of the Spacetime Scalar Curvature 

The Einstein equation contains the spacetime scalar curvature tensor ^R so we are inevitably forced to 
find its 3+1 decomposition. Otherwise, the 3+1 decompositions of the Einstein equation can not to be 



constructed. Therefore, let us take the trace of the equation ( |147[ ) with the dual induced 3-metric y"^: 

1 



r"^/.//^ %v = -^y^£mKij - ^DiD^N + R + K^- IKijK'i , 



(149) 
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it reduces to 



r 



N N 



(150) 



The 1*' term of the equation (150 1 contains y^'' and R^iv In order to take the trace of R^y, we need to 



rewrite y'^" in terms of which is given in the equation ( 35 1 



(151) 



Observe that ( 150 1 is a scalar equation. Because of this, we need to find the explicit form of the 2"'^ term of 
( [Bol l: Now, 



(152) 



The last term of the equation ( 152 1 contains the Lie derivative of the dual induced 3 -metric, 

-Lmj'K along the normal evolution vector m but we do not know what JLiaj'K directly. Fortunately, we 



know the Lie derivative of the induced 3-metric along m [equation (71 1]. Therefore, to find XmT'^> we will 
use JLmYij by starting from the relation jikj'^^ - : 



Let us multiply the equation ( 153 1 by y'' 



I ^ i,i ^ k,k ^ I 



(153) 



= iNy^y^Ku 



With this result, the equation ( 152 1 turns into the form that we want 



r'^Xm^/y = LmK - INKijK'j 



Finally, by substituting the results of (151 1 and ( 155 l into the main equation (150 1, we get 



'^R + "^Ruvn^n" = --£mK - -DiD'N + R + K^ . 



N 



N 



(154) 



(155) 



(156) 
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To find the 3+1 decomposition of the spacetime scalar curvature we need to rewrite the term in the 



equation ( 156l which contains the spacetime Ricci tensor. This is done with the help of the scalar Gauss 



relation ( 128 1: 



%,ti'n' = + -K' - -K,iK" ^ 



Thus, by substituting the results of ( 157 1 into the equation ( 156 1, we will get the 3+1 decomposition of 
spacetime Scalar curvature as 

^R = R + K^ + KijK'^ - ^£rnK - ^DtD^N . (158) 
N N 

IV. THE 3+1 DECOMPOSITION OF EINSTEIN EQUATION 

In this chapter, we will suppose that the spacetime that we are going to deal with is a solution of the 
Einstein equation with zero cosmological constant. Now, the Einstein equation: 

^R-^^Rg^SnT. (159) 

It can be assumed as an equation of machines where each machine has two slots for inputs. Alternatively, it 
can be written in terms of matter parts 

4R = 87r[T-^rg]. (160) 

where T stands for the second rank stress-energy tensor and T is its trace. Then, by substituting the corre- 
sponding inputs into the slots, one can reach the 3+1 decomposition of the Einstein equation. In order to 



find it, as we see from the equations ( |159| ) and ( |160| ), we need first to find what the 3+1 decomposition of 
the stress-energy tensor T: 

A. The 3+1 Decomposition of the Stress-Energy Tensor 

Now, the full projection of the T along the unit normal vector " is nothing but the energy density E 
evaluated by the observer, 

E = T{n,n). (161) 
And the full projection of T onto is bilinear form stress tensor S: V(u, v) € TpiJ^i), 

S = /t = T(u,v), (162) 
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or in components form 

Sap^y^aJ^pT^jtv (163) 



What the equation ( 163 1 says is that the action of the component of the force along Ca onto the infinitesimal 
surface of the normal along ep is given by S aii- Here, Ca and Ca are the spacelike vectors with respect to the 
frame which does the measurements. 

Finally, the mixed projection of T is the component of the 1-form momentum density p: That's, 



p = -T(n,v), (164) 



or in terms of the components 



Pa = -n'^y^'aT^v (165) 



Now, let us think reversely. We naturally expect that the stress-energy tensor must be the union of the 
corresponding projections of it: That's, 

T = Sn^n + S -i-p(8)n + n(8)p. (166) 



From the equation ( 166 1, it is easy to show that the important relation between the traces of T (i.e. T) , S 



(i.e. S) {with respect to g] and E is given by 

T = S-E. (167) 

B. The Projection of the Einstein Equation 

With the help of the fundamental relations of the 3+1 formalism that we have found in the previous 
chapter and the 3+1 decomposition of the stress-energy tensor, we will find the 3+1 decomposition of 
the Einstein equation which are known as the dynamical Einstein equation, the Hamiltonian and the 
Momentum constraint equations. 

1. Full Projection onto 



First of all, let us start by fuUy projecting the Einstein equation of ( |160[ ) onto the hypersurface If! 

^*4R = 87r(-^*T-irt*g). (168) 
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Now, the knowledge of T is equal to the difference between S and E (see the equation 167 1, the full 



projection of the T onto E, (the equation 162 1 is S and the full projection of the g is the induced 3-metric y 



with the equation ( 148 1 of 



^ --XmK - ^DDA^ + R + /:K - 2K K 



leads us 



- ;^^mK - -^DDA^ + R + - 2KK = 87r(s - ^[5 - E^y) , 
i-X^K = --^DDA^ + K + KK- 2K.K + %n(^[S - E]y - s) , 
XmK - -DDA^ + NR + KK- 2K.K + 4n{[S - E]y - 2S) , 
or in the components form 

-CmKa)3 = -DapDapN + N\Ral} + KKa/3 - IKa^K^ p 

+ An{[S -E]y,p-2Sap)]. 



(169) 



(170) 



Since the equation ( 170 1 is totally composed of the quantities of E^, we can rewrite it in terms of the spatial 
indices, 



Lry^Kij = -DiDjN + N[Rij + KKtj - IKn^K^j 

+ 47:{[S -E]yij-2Sij)]. 



(171) 



The equation ( 171 1 is known as the dynamical part of the Einstein equation and it has 6 independent 
components. 



2. Full Projection Perpendicular to 



Secondly, let us fully project the Einstein equation ( 159 1 along the normal unit vector n 

'^R(n,n) - ^'^Rg{n,n) - 87rT(n,n) . 



From the equation ( 161 1, the component form of the equation ( 172 1 becomes 



(172) 



(173) 



By using the scalar Gauss relation ( 128 1, we get the Hamiltonian constraint part of the Einstein equation 



as 



R + K^-K,iK'J = l6nE. 



(174) 



Since the Hamiltonian constraint equation {174) is a scalar equation, it has only 1 independent component. 
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3. Mixed Projection 



Let us project the Einstein equation ( 159 1 either one-times onto the hypersurface or one-times orthogonal 



to It 



4R(n,^(.)) - ^-'Rg(n,^{.)) - 87rT(n,t(.)) ^ 'R(A, r(.)) - 87rT(n,^(.)) . 



From the equation ( 176 1 we know that the mixed projection of the T is p. By using this knowledge, the 
component form of the previous equation becomes 



(175) 



Now, by substituting the Contracted Codazzi Relation ( 136 1 into the equation ( 175 1, we get the momentum 
constraint part of the Einstein equation (or simply Momentum constraint equation) as 



(176) 



where ( 176 1 has 3 independent components. Moreover, it can also be rewritten in terms of the spatial indices 



DjKJi - DiK = ^npi . 



{Ill) 



C. The 3+1 Dimensional Einstein Equation as a PDEs System 



Now, the coordinate systems adapted to the flow of the hypersurfaces (96 1 allow us to transform the 3-1-1 
Einstein equation into a set of partial differential equations (PDEs). 



1. Lie Derivatives Along "m " as Partial Derivatives 



1. For£mK: 



The equation (98 1 provides us to decompose the Lie derivative of the extrinsic curvature along m in 
terms of the Lie derivatives of K along dt and along the shift vector yS 



XmK — X.^ K — X^K . 

d, P 



(178) 



The equation ( 178 1 says that K is also a tensor field of the hypersurface. We wish to write the 



tensor components relative to the well-defined coordinates. Then, the K turns into the partial 
derivate relative to the time coordinate t 



Or 



dt 



(179) 
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Similarly, the Lie derivative of K along jS turns into the partial derivatives relative to the spatial 
coordinates 



H " ^ '^'^ ''^ ■ 



5x' dx3 



(180) 



2. For£mr: 



Identically, the equation (98 1 provides us to decompose the Lie derivative -CmT as 



(181) 



Again, because of the same reason, the Lie derivative of the induced 3-metric becomes the partial 
derivative. 



dy. 



(182) 



and also 



= ykjDi/3'' + yikDj/3'' 
= Di/3j + Dji3i. 



(183) 



With the help of the form of the corresponding Lie derivatives ( 178 1 and ( 181 1 in the foliation adapted 



coordinates, we get the 3+1 decomposition of Einstein equation within these coordinates 

d 



(184) 



- - £^ ]Kij = -DiDjN + N\Rij + KK^j - IKi^K^j + 4n[{S - E)yij - ISij] \ , 



(185) 



R + K^-KijK'j = l6nE, 



(186) 



DjKii - DiK - %npi . (187) 
This set of equations is known as the 3+1 dimensional Einstein system. 
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V. CONFORMAL DECOMPOSITION OF THE 3+1 EINSTEIN EQUATION 

In the Ricci and Cotton flows Il23l . EH, the Riemannian metric is being mapped by particular difi'eo- 
morphisms. Furthermore, if the mapped metric is a scale times the initial metric then it is called the gradient 
Ricci and the gradient Cotton solitons. As in the Ricci and Cotton flows, we will examine the evolution of 
the hypersurfaces as if there is a well-defined Riemannian conformal background metric y and it is smoothly 
mapped into the 3-metric of Z, by a positive scalar field *F during the flow. This type of the flow of the metric 
is known as the conformal transformation (or flow) of the induced metric of ( see Lichnerowicz [61) 



r-T^y. (188) 



Now, let us take a look at what Eric Gourgoulhon [Ij says about the York's publications |[T4l. ifTSl which 
give the importance of the conformal transformation in gravity : " In 1971 - 72, York has shown that 
conformal decompositions are also important for the time evolution problem, by demonstrating that two 
degrees of freedom of the gravitational field are carried by the conformal equivalence classes of the 3- 
metrics. A conformal equivalence is defined as the set of all metrics that can be related to a given metric jij 



by a conformal transformation like the relation of ( 188 1". Now, we know that the Weyl tensor gives whether 



a given spacetime is conformally flat or not [that's, the background metric in the equation (188) is flat] and 
it is valid for the spacetimes whose dimension is greater than 3. Here, in the lower dimensional cases the 
conformally invariant Cotton tensor, Cijk, lfT3]l takes the role of the Weyl tensor: 

Cijk - Dk{Rij - ^Rjij) - Dj(Rik - ^-Rjik) . ( 1 89) 

Moreover, the York ifT?! . |[T5]| showed that the C'J is another conformally invariant tensor, 

dI = ^I^C'i, (190) 

where the C'^ is the well-known Cotton- York tensor |[T4]| . |[T5l . |[T3]| which is constructed from the Cotton 
tensor pSOll US as 



C'7 . _ 1 '-^/c,,,,,/"^' = e"^'Dk(R^i - \R5i) , (191) 



which satisfies the following properties 
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1. Symmetric: 



(192) 



= 0. 



2. The traceless: The Cotton- York tensor is traceless, then, let us first rewrite this in terms of the sum 
of the Einstein tensor: 

Ci = L''<'D,(Rji - jRdJi) + y"D,(R', - jR6'i) 



4^" 



1 



-^e'''D,&i + y''DtG': 



By taking the trace of ( 193 1 with respect to yij: 

1 



y.jC'J = -e'"DkGu + -e}"DkGj,. 



Because of the symmetric and antisymmetric relations in the equation ( 194 1, we get 



(193) 



(194) 



yijdj = 0. 



(195) 



3. Divergence-free, (i.e. transverse): DjC'J - 0. 
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Now, let us prove that the tensor defined in the equation {190 \ is really conformally invariant under the 



(196) 



transformation {188) by using the following particular conformal transformations: 

2. detijij) ^ {^>ydet{yij) ^y = ^'^y . 

3. e ^ ^ e'^' = vF-^e'*^' where y - det{yij) . 

4- Cynkl ~ Cmkl ■ 

Then 

From here, we are going to find what happens to the 3+1 form of the Einstein equation under the 
conformal transformation. Therefore, in order to find the explicit form of the conformal 3+1 expression of 
the Einstein equation, we first need to find the conformal form of the fundamental quantities: 

A. The Conformal Form of the Intrinsic Quantities 

Because the determinant of y depends upon the choice of coordinates so the conformal factor is not a 
scalar field. Thanks to the flat background metric, we achieve to make the conformal factor a scalar field 

m. 

1. Conformal Connection 

Now, suppose that there is a conformal background metric y on the hypersurface whose intrinsic con- 
nection is metric-compatible 

by^O. 
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And associated ChristofFel symbols are defined as 



Here notice that the partial derivatives are taken with respect to the coordinates (x'). Since the intrinsic 
covariant derivative of any tensor field T can be taken relative to two connections D and D, then, it is 
normal to expect that there must be a transformation between them. And this transformation is given by 



. (198) 



<- kjr^ Jl-l-Jq ' 



Here C'^ij - - ft. Also, it a tensor field because it is defined as the difference between the Christoffel 
symbols. And its explicit form is 



/^k _ rk _ f^k 
^ ij - ^ ij A ij 



■k pm 
ij " mi ij 



- [ - Dmj - rUrmj - tTjyn.] - 2tp^i] 



(199) 



Because of the conformal transformation of induced 3-metric and its dual. 



ro = *%-, (200) 



,ij - vp-4^i7 



(201) 
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We can rewrite the equation ( 199| ) fully in terms of the conformal quantities as 



C^7 = l7'''{Dmj + Dj7il - blJij 

= 2{(5^yA(lnY) + 5^■Dy(ln'^') - D*^(ln'F)y;y 



(202) 



The equation ( 202 1 is playing an important role in the conformal transformation. As a sample, let us find 



the relation between the v € r(S,) 



Djv' = Djv' + C'jkv^ 



= bjv' + 2^^6'jbk{\n ^) + v'Dy(ln - v^yjkb\\n ^)| . (203) 



Let us do the change of j — > / in the equation (203 1, 



Div' = biv' + 2|v*(5',D;t(ln^) + v'D,(ln4') - v*=ri^D'(ln »P)| 
- Av' + 2(3v^D;t(ln*F) + v'D,(ln^) - v''bk{ln^')] . 



k ^ i k ^ i 

Thus, we get the conformal form of the divergence of the vector tangent to E as 

Div' - Av' + 6v'A ln4' 
= '¥-^biC¥^v') . 

2. Conformal Transformation of the Intrinsic Ricci Tensor 

1 . Conformal Relation of the Ricci Tensors in terms of the Tensor Field C 
The corresponding Ricci identity is 



(204) 



{DiDj - DjDi)v^ = R'^iijv' . (205) 



By doing the suitable operation of contraction and change of dummy indices, the equation ( 205 1 
becomes 

Rijyj - DjDiV^ - DiDjv^ . (206) 
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Now, with the help of the general transformation relation ( 198 1, the equation (206 1 can be written in 
terms of the conformal tensors of R , C and the conformal covariant derivative of C, 

Rijvj - DjiDivJ) - C'^jiDkvj + djuDiv'' - Di{Djvj) 

- Di\Djvi + djkv''] 

- DjDivJ - DiDjvJ + bjdikv^ - djidkiv' + Cjjkdiiv' - Did jkv^ 
= Rijyj + Djdikv'' - djiduv' + djtduv' - bidjkv'' 

j 1 ^ j j k^ j 

= RijvJ + Dtdijyj - dudkjv^ + dikdijyJ - Didujv^ . 

Because is an arbitrary vector field, we get 



(207) 



■^^'7 = -^'7 + ^kdij - bidkj + dijdik - dudkj ■ 



(208) 



2. The Conformal Transformation of the Intrinsic Ricci Tensors in terms of the Conformal Factor 
The conformal equation ( |208| ) can be rewritten in terms of the Conformal factor ^ by using the 



equation of (202 1. For simplicity, let us work on the terms of the equation (208 1 which contain the 
conformal covariant derivative of tensor field C : 



(a) VoxDidkj 



dij = 2 (5^■D/ln'^') + S'^jDiiln^) - D\ln^)yij 



Let us do the change / ^ ^ in the equation ( 209 1, then, it becomes 

dkj = 2{(5\DXln^) + 5^DKln^) - D\ln^')%j 
= 2|3D/ln'F) + D/lnf) - D/ln^) 
= 6D/ln^), 

so we get 



(209) 



(210) 



(211) 
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(b) For buC\j : 



= 2 D/D/ln'F) + DjDi{\nm) - DkD\\n'¥)yi 



(212) 



Let us substitute the results of (211 1, (212i and the explicit formula (202 1 for the components of C 



into the main equation (208) 



= Rij + 4D,D/lnT) - 2D^D*(ln^)7,y - eD/D/ln'F) 



+ 2 (5%D/ln »P) + d^PiiXn ^) - D*^(ln x 6Dk In 



- 4|(5^•D,(ln + 5^D,(ln 4^) - D'^(ln 4')y,7|x 
{^'iD/ln + 5',D^(ln - D'(ln 4^)^; 
= - 2D,D/lnT) - 2D^D*^(ln^)f,-y 



- n{& In ^)(D^ In »P)7;y - 43^6' k0i In T)(Dy In 

- 45^•5'y(A In ^'){Dk In »P) + 4(5^•(D; In 4^)0' In 'F)riy 

- Ad''i6\{bi In T)(Dy In »P) - A5''i6^j{bi In »P)(D^ In ^F) 
+ 4(5^0,- In ^)(D' In + 4(5'^0'' In In T)y,7 
+ 4(5'/D^ In ^)(D^ In »P)y,7 - 4(0^^ In Y)(D' In , 



(213) 



by collecting the identical terms of the equation (213 1 in the each corresponding clusters, we get 
conformal transformation of the Ricci tensor as 



Rij = Rij - 2D,Dy(ln'F) - 2D^D'^(ln ^)yiy + (A ln'F)(Dyln4') 



(214) 



-4(D;tlnT)(D'^ ln'F)y, 
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3. Conformal Transformation of the Scalar Intrinsic Curvature 



Let us first take the trace of the equation (214 1 with respect to the dual induced 3-metric y'^: 

R = yjRij 



= . - IfJOiDjiln - 2y'^yijDkD\ln ^) (215) 

+ 4fj{Di In 'F)(Dy In ^) - 4f^fij{Dk In ^)(D*^ In T)] 
= Y""^!^ - 8[d,D' In + (A In ^)(D' In Y)]] 

Here we need to modify the term DjD' In *F of the equation ( 215| : 



- ^'-^DiD'^' - {Di In 'F)(D' In . 



(216) 



Thus, the substitution of the equation (216 1 into the equation (215 I results in the conformal transformation 
of the intrinsic scalar curvature of 

R ^ - 8^"^ AD'*i' • (217) 

4. Conformal Transformation of the Extrinsic Curvature 

Since the trace and traceless parts of the 3+1 Dynamical Einstein equation transform differently un- 
der conformal transformation, we need first to decompose the extrinsic curvature into the trace part and 
traceless part. 

1 . The Extrinsic Curvature in terms of Trace and Traceless Parts 
Now, the traceless part of the extrinsic curvature is defined as 

A-K-^^y, (218) 

such that fryA = y'M;y = 0. Therefore, the covariant and conravariant components of the extrinsic 
curvature can be rewritten in terms of the trace and traceless parts. 



Kij = Aij + ^Kjij and K'j = A'^ + ^Ky'j . (219) 



2. Conformal Transformation of the Traceless Part 
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As in the conformal transformation of the induced 3-metric jij, the conformal transformation of the 
traceless part of the extrinsic curvature must be something like 



(220) 



We will see that the choice a = -4 gives us the evolution equations for the conformal factor ^, the 
conformal 3-metric jij and its dual. On the other hand, the choice a = - 10 gives the conformal form 
of the moment constraint equation. 

(a) For the 1*' choice of a = -4 



The Lie derivative of the induced 3-metric is given by the equation (71 1. Let us find what 



happens to it under the conformal transformation by using the equations (200 1 and (219 1: 

ZmYij = -2NKij 



(221) 



then 



"V^-Cmfij + (Xm*P^)r,-; = -2NAij - -NKjij , 



- -INAij - ^NfC¥%j - (Xm^')r,7 
£m7ij - -2N'¥-%j - \NKyij - ^^(inT^J-.. 



(222) 



-2N'i>-^Aij - -NKyij - 4Xm(ln»P) 



7ij ■ 



Therefore, the equation ( 222 1 becomes 

£m7ij - -2N'¥-%j -1{nK + 6£m{ln'¥)\yij . 



(223) 



Since the Ajj is traceless, let us multiply the equation (223b by the conformal dual 3-metric y'^: 



f'-Cm7u = -2NW^f^Aij - 3{a^^ + 6i;m(ln>I')}r%- 
= ~iNK + 6£„,{ln^')\x3, 



(224) 



so we get 

y'-Cm7ij = -2[nK + 6£™(ln»I')} . 
Now, the variation of the determinant of an invertible matrix is given by 

6[lndet:n) = tr{^~^ x 6J{) . 



(225) 



(226) 
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Let us do the changes of J?l ^ y,y and 6 Xm in (226 1 so it turns into the left hand side of the 



equation (225 1. Then, the left hand side of the equation (225 1 becomes the Lie derivative of 
a scalar field along m which allows us to decompose the Lie derivative along m into the time 
derivative and the Lie derivative along the shift vector /3 



d 



(227) 



Because the time derivative of the scalar field In detifij) vanishes, the equation ( 227 1 reduces 
to 

£mln detifij) ^ -Zp\ndet{yij) = - f^Zpjij 

= -fj{/3''Dkyij + nyD,/ + yikDj/3''] 

= -r{ykjDi/3' + yikDj/3'} (228) 

= -S'kDifi'' - Shop'' 

= -2D0 . 



Furthermore, by substituting the result of (228 1 into the equation (225 1, we get the evolution 
equation for *F under conformal transformation as 

(--£p)\n^= -{Di/3' - NK) . (229) 

By inserting the equation ( 229| ) into the equation ([223 1, it becomes 

(230) 



2 ~ 



For consistency in the equation (230 1, we must have such a conformal transformation of the 
traceless part of Kjj as 



Aij - *P Aij , 



(231) 



which says that the corresponding a must be -4. Therefore, we find the evolution equation for 
conformal metric during the conformal transformation as 



d 



(j--Cls)7ij--2NAij--D,/3'yij 



(232) 



"■dt 3' 

Now, the conformal transformation for the contravariant component of the traceless part of 
A'^, is obtained by 

pi] = 

= ^-^A'j . 



(233) 
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Finally, let us see how the dual conformal 3-metric y'j evolves under the conformal transfor- 



mation by starting from the equation ( 232 1 : 



Thus, we get 



(j-£/^)f' = 2NA'J + lDtP'y'J 



(234) 



(235) 



^dt 3 
(b) For the 2"^ Choice of a - -10 

In this case, we will start with the decomposed form of the conravariant extrinsic curvature K^^ 



[see the equation (219)] and take the divergence of it. That's, 

^'V = A'^' + ^-Ki '^ ^ Bit' - D/A'^' + -D'/T. 
Z ' ' ^3 



(236) 



Now, the equations of ( 198 1, (202 ) and (21 1 1 provide to rewrite the term D^A'^ of the equation 



(236 1 in terms of the conformal quantities, 

DyA'-'' - DyA'-'' + CjkA^'^ + d jkJ^^ 

- DjA'i + l\5' jbk In *F + 5\bj In »P - D' In "^jj^A^^ 
+ 6Di:ln»PA* 

= DyA'^' + IS jA^'bi, In ^ + 15\A^'bj In ^ - IjjkA^'b' In 
+ 6Di In^A'*^ 

- bjA'' + 2A'*^D^ In Y + 2A'^Dy In *F - T^'^'y j],A^^& In *F 

+ 6Di ln*FA'^ 



(237) 



Because A*^^ is traceless, the related term vanishes. Then, ( 237 1 turns into 

DyA'^' - DyA'^' + 2A'*^Di. In T + 2A'^Dy In + 6D^ In ^PA''^ , 

j j 

with the given change of indices, we get 

DyA'^' = bjA'j + lOA'^Dyln^ 

- xp-10Dy(Y^°A'^) . 



(238) 



(239) 
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Notice that for consistency of 



(240) 



the corresponding a must be -10. 



Finally, let us insert the equation ( 343 1 into the momentum constraint equation {177 \ 



.. 1 . .. 2 ■ 

DjA'^ + -D'K - D'K = %np' ^ DjA'^ - -D'K - ?,7Tp' 



And let us insert the conformal form of the DjA'J (239 1 into the equation (241 1 



Now, the modification of 

DjK = DjK ^ y^'DjK = D'K - y'JDjK 



provides us to rewrite the equation ( 242 1 as 



Thus, the the conformal transformation of the momentum constraint is 

Dj[A'j) - ^"V^D'K = S/r^iV' • 

Finally, let us find the conformal transformation of Afj 

A,, = y'yJ'A'' 

= {^'-\''){^'-'yJ'p''A") 
= ^y^yj'A'"' , 



then, we get 



(241) 



(242) 



(243) 



(244) 



(245) 



(246) 



Aij - ^ Ajj 



(247) 



B. The Conformal Form of the 3+1 Dimensional Einstein System 



As we said before, the trace and traceless part of the 3+1 dynamical Einstein equation ( 171 1 transform 
separately under conformal transformation. Therefore, we need to first deduce the related decomposition of 
it. And then, we will be ready to construct their conformal forms. 
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1. Trace and Traceless Parts of the 3+1 Dynamical Einstein Equation 



1. Trace Part of the 3+1 Dynamical Einstein Equation 



Let us take the trace of the 3+1 dynamical Einstein equation ( 171 1 with respect to y 



y^£mKij = -DiD'N + nIr + K^- IKijK'j + 4n[s - 3E) 



From the equations ( 155 1 and (248 1, we have 



£mK = yj£„,Kij + INKijK'j 

= -DiD'N + nI^R + K^- IKijK'^ + 47r(5 - 3£)| + INKjjK'j 

= -DiD'N + nIr + K^+ 47t[s - 3e)\ . 



Now, the Hamiltonian constraint equation ( 174i 



R + K^- KijK'j = \6nE ^ R + = \6nE + KijK'j , 



(248) 



(249) 



(250) 



by substituting the previous modification into ( 249 1, we obtain the trace part of the dynamical 3+1 
Einstein equation as 



£^K = -DiD'N + NiKijK'J + 47t[s + e) 



(251) 



2. Traceless Part of the 3+1 Dynamical Einstein Equation 



Now, let us now decompose the left hand side of the 3+1 dynamical Einstein equation ( 171| ) by using 



the equation ( 2 



-CmKij - £m{Aij + l^Kjij) 

= £mAij + ^{£niK)yij + l^K[£^yij) 
- £mAij + ^{£mK)yij - ^KKij , 

where we used (71 1. Therefore, the Lie derivative of the traceless part Aij along m is 

£mAij = £mKij - U£mK)yij + ^KKij . 



(252) 



(253) 



Notice that the first term on the right hand side of (253 1 is the 3+1 dynamical Einstein equation ( 171 



61 



and the second term is the trace part of it ( 249 1. Then, the explicit form of ( 253 1 is 
Lr^Aij = -DiDjN 

+ N^Rij + KKij - IKiuK^j + 4n[{S - E)yij - 25 

-U- DkD^N + n[r + K^ + 4n(s - 3E)\\yij 



2N 



- -DiDjN 



(254) 



+ N{Rij + \KKij - IKi.K'^j - \K^yij - ^n(s ij - Uyii) 



+ -\DkD'N-NR\yij. 



Since we wish a totally traceless equation, we need to get rid of the terms that contain the K : 
^-KKij - 2Ki,K^j - ^-K^yij = ^-K{Aij + ^-Kyij) 

-2(Ai, + lKyi,)(A'j + \K6'^j) 



= -KAij-2AikA'j. 



(255) 



Thus, by substituting (255 1 into (254 1, we obtain the traceless part of the dynamical 3+1 Einstein 
equation as 

£„,Aij - -DiDjN + Ni^Rij + ^KAij - 2AiuA'' j - u(s ij - \syii)] 
+ ^\D,,D^N-NR^yij. 

After constructing the conformal transformation of the fundamental tools, it is time to construct the confor- 
mal transformation of the 3+1 Einstein equation : 



(256) 



2. Conformal Decomposition of the Trace and Traceless Parts of the Dynamical 3+1 Einstein Equation 



We are trying to construct the corresponding time evolution equations. Therefore, we need use the 
a = -4 case. 

1. Conformal Form of the Trace Part of the 3+1 Dynamical Einstein Equation 



The trace part of the 3+1 dynamical Einstein equation (251 



IS 



d 



(- - £f!)K - -DiD'N + N\KijK'j + 47r[s + e) 



(257) 



For simplicity, let us find only the conformal form of the terms DiD'N and KijK'-i separately: 
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(a) The conformal form of the term DjD'N 



We have found that the conformal transformation of the divergence of a vector v ( 204 1 is given 
by 



Since the gradient of a scalar field is a vector field, we will take v' = D'N, 



(258) 



v' = D'N = y'^DjN = 'f-'^yJDjN = '¥~''D'N . 



-4 iSi , 



(259) 



Let us substitute (259 1 into (258 1, 



DiD'N = 'i'"^Di{^'^D'N) 

= "^-^bii^^bi&N + 2TO,'FDW) 
= "V-^ipi&N + Ibi In tow) . 



(b) The explicit decomposition of the term KijK'J 



Let us use decomposed form of the extrinsic curvature (219 1 



Because the trace of A,-, is zero, (261 1 reduces to 



= AijA'j + -K^, 



(260) 



(261) 



(262) 



where we used the conformal transformations (231 1 and (233 1. 



Finally, by substituting ( 260| ) and (262 1 into the fundamental equation (257 1, we will get the confor- 
mal form of the trace part of the 3+1 dynamical Einstein equation as 



(— - £/j)K = -'¥"^[bib'N + 2biln^'&N) 

+ N[AijA'j + ^K^ + 4n[E + s)] 



(263) 
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2. Conformal Form of the Traceless Part of the 3+1 Dynamical Einstein Equation 



The traceless dynamical 3+1 Einstein equation (256 1 is 



+ -I^DtD'N-NR\yij. 



First, let us find the conformal form of the terms Xm^/y and DiDjN 



(a) The conformal form of 



(264) 



= 'P^XmA/y + 4^'^{£n.'¥)Aij 

= 'i"'£rnAij + 4^'\£^ln^)Aij 

= Y^XmA/y + 4»P4(1[D,/ - NK])a 

= 'I'^jXmA/y + y^Dkli'' - NK)Ai}^ , 



(265) 



where we used the conformal evolution equation of ^ {229 1 and the conformal transformations 



{231 1 and {233 ). Finally, let us use the equation (256 1 of XmA,y in the following equation 



XmA,y - "y-^ZmAij - -(d^/ - NK)Ai 



+ 



N[Rij + ^-KAij - 2Ai,A'j - %n[Sij - ^^y^y]) 



+ -[DkD''N-NR^yij 



-:(Dk/3^ - NK)Aij . 



(266) 



(b) The conformal form ofDjDjN 



DiDjN = DiDjN 



= DiDjN - C^ijbkN 



(267) 

DiDjN - 2[biNbj In T + DjNDi In T - bkN& In "^jij] . 



= bibjN - 2[5^■Dy In + d'^jbi In x buN 
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Finally, let us substitute the equations (214 1, (260 1, (266 1 and (267 1 into the equation of (264 1 
Z^Aij = Y"'^! - iOiDjN - l[biNbj In Y + DjNDi hi 'F - DkNb'' hi ^y,y) j 



+ N\\Rij - 2D,Dy(ln^) - 2D^D''(ln ^)y, 



+ 4(Di In »P)(Dy In ^F) - A{Dk In '^){& In 'V)yij 
+ ^-K^y'Aij] - m'^fAikAj, - 87T[Sij - ^S^'^jfi (268) 
+ ^i^^'-'^^Dk&N + 2bk In "HD^n) 

- N'V'^R - %[bk& In ^ + {bk In In 'F)])|'FV^j 
-^[bul3''-NK\Aij. 

After some algebra, we reach the conformal transformation of the traceless part of the 3+1 dynam- 
ical Einstein equation as 

-CmAij = -^bk/3%j + Ni^KAij - 2f%kAji - S7r[¥-^Sij - ^57,7]} 

+ ^P""*/ - bibjN + ibiNbj In Y + ibjNbi in »P 



+ 



+ 



^ [DytD^W - 4bk In TO'^A^Jr/y 



2bibj In 4^ + 4D, In TOy In 



'-\bk& In 4' - 2Di In In 'Fjy^y] 



(269) 



3. The Conformal Transformation of the Hamiltonian Constraint 



The Hamiltonian Constraint equation ( 174 1 is 

R + K^-K,iK'j = \6jtE. 



(270) 



Now, let us substitute the conformal transformation of R (217 1 and the decomposed form of KijK'J (262 1 



into the Hamiltonian Constraint equation (270 1, 



^P""^^ - iW^b0'¥ + K^- AijA'j - -K^ = l6nE . 



From the equation (27 1 ), we get the conformal transformation of the Hamiltonian constraint as 

bi&'V - -R'V + l-AifA'j - —K^ + InEW^ = . 



Due to the relation of AijA'^ = Y ^^AjjA'-', the equation (|272| becomes 



8 8^ I 12 ' 



(271) 



(272) 



(273) 
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which is known as Lichnerowicz Equation 161 Q. 

Thus, the conformal transformation of the 3+1 dimensional Einstein system ^ can be summarized 

as, 

(| - X/j) In * - ^^(Dfi' - NK) , (274) 
(- - £is)7ij - -2M,-y - -D,/3'y!j , (275) 



(— - £/j)K - -^■'^(d,DW + 2D, In tow) 



(276) 



+ ^P-^^I - D;DyA^ + IDiNDj In 4' + 2DyA^Di In Y 

+ ^ [D^tD^A^ - Abk In TO*A^]r,y (277) 

+ N[Rij - ^Rjij - Ibjbj In ^ + Abi In TOy In Y 

2 1 
+ ^(^-^^^ In »P - 2D;t In TO*^ In 'P)?/;]} , 

D,D'^ - ^CT + ^A,yA'^*F"^ + [27r£ - ^K^] = , (278) 
bj[A'^) - ^"^^b'K = Stt^^V ■ (279) 

C. The Isenberg- Wilson-Mathews Approach to General Relativity(IWM) 

In IWM model lITTI . lITS! . the spacetime is assumed to be foliated by a continuous set of (Lt)te'R such 
that the foliation is maximally sliced (K = 0). Here, the induced 3-metric is conformally flat which means 
that its conformal background metric is flat, 

yij = fij- (280) 
(280 1 is implies that the Cotton- York tensor |[T4ll . |[T5]| . lfT3l vanishes. Furthermore, the conformal Ricci 



tensor is zero. Thus, the conformal 3+1 Einstein equation turns into, 

(|-^/j)lnY=^A^', (281) 
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- -^'-'^{DiD'N + 2Di In "VD'n) + N[AijA'j + 47r(£' + 5)) . 



(282) 
(283) 



d '2, ( 1 



+ vj/-4| _ DiDjN + IDiN Dj In ^ + IDjN Di In Y 
+ ^[DkD'^N - 4Dk In^" D''N]fij 
+ n[- IDiDj In *P + ADi In In »P 



+ 



|(DiD^ In 'F - 2D^ In WD^ In , 



(284) 



(285) 



DjA'^ + GA'^Dj In T = 87r^'*p' . 



(286) 



Here the equation (286 1 is obtained by using the relation AjjA'j - Y ^^AjjA'^ in the momentum constraint 



equation (245 I. In order to find the IWM conformal system, let us work on the equation (282 ) : Because of 



the metric-compatibility, we have 



= fkjDi0' + fikDjpK 



Now, since time derivative of the fij is zero, the equation ( 282 1 turns into. 



Let us multiply (288 1 by /™/- 



(287) 



(288) 



By change of the indices m ^ i, n ^ j, we can rewrite the equation (289 1 as 

1 



A'J 



IN 



where 



(289) 



(290) 



(Ipf = D'|i^ + D^|3'--Du|3'f~ 



{29\) 
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is known as the the conformal Killing derivative operator. Moreover, with the help of ( 290 1 the corre 



sponding momentum constraint equation (286 1 can be rewritten as 



A + l^D'Dj/3j + 2A'j[6NDj In ^ - DjN) = IGnN^^p' . (292) 
Thus, we get the conformal IWM system as the set of 

A + 2Di In ^'D'N - n[47:{E + S) + AijA'j)^^ , (293) 

A »P + [ ^AijA'j + InE]^'^ - , (294) 

A p' + \^D'Djl3^ + 2A'}[6NDj In - DjN) = \6nN'i''^p' . (295) 

VI. ASYMPTOTIC FLATNESS AND THE ADM FORMALISM FOR GENERAL RELATIVITY 

A. The Asymptotic Flatness 

In this chapter, we will deduce the conserved quantities of the ADM mass, linear momentum and an- 
gular momentum of a given hypersurface E^. Since these quantities can be only in the globally-hyperbolic 
asymptotically fiat spacetimes (i.e. the spacetimes which approaches asymptotically to the well-defined 
spacetimes such as the Minkowski, AdS). Therefore, let us first see review what the asymptotic flatness is: 
The asymptotic fiat spacetime is such a particular spacetime for the massive objects in which it is assumed 
that there is nothing in the universe except these objects. Now, a globally-hyperbolic spacetime is called 
asymptotically flat if each of its Cauchy surface has a background metric f with signature (+,+,+) such that 
f is fiat, can be diagonalized in a particular coordinate system on the HI, |[T4l . Moreover, in the case of 
spatial infinity, r ^ oo, the decay of jij and their spatial partial derivatives must be something like 

jij = fij + 0[r~'] , (296) 

= OVr-^^ . (297) 
And also as r — > cx3, the decay of Kjj and their spatial partial derivatives must obey 

Kij = 0[r-^] , (298) 
dKu , 

= 0[r'^^ ■ (299) 
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B. The Hamiltonian Formalism for the General Relativity 

The Hamiltonian model approaches a physical state at a certain time and gives the evolution of the state 
as t varies. This model is being transformed into the gravitational theory as a state on a particular spacelike 
hypersurface. Now, the gravitational theory is a covariant theory and locally has Lorentz symmetry. The first 
attempts tried to start with the spaceUke hypersurface that is free of choosing coordinates to avoid breaking 
of the crucial properties of the gravitational theory ||9l, ifTOl . However, it is then hard to define initial state of 
practical problems. In order to write the Einstein equations into the Hamiltonian form, people started to give 
up the main properties of the gravitational theory by choosing a family of particular coordinate systems such 
that " x" = constant" corresponds a spacelike hypersurface. Contrary to the unknowns {jij , Kij , N , fi' ) in 
the PDEs form of 3+1 Einstein system, Amowitt, Deser and Misner have proposed the ADM formahsm of 
the General Relativity II16I1 in which conjugate momentum of the induced three-metric jij, n'J - -s/yiKy'^ - 
K'j), is used instead of Kjj. Moreover, in the ADM formalism, n'^ and yij are the dynamical variables and 
the Lapse function N and the shift vector yS are taken as Lagrange multipliers |[T6l . 

In this section, we will first deduce the corresponding Hamiltonian form of the vacuum field equation by 
mean of the 3+1 decomposition of the spacetime metric that we have found in the V chapter (the equation 



1 15 1 and the knowledge that the boundary term is zero. Secondly, we will deal with the general case. That's, 
we will deduce the corresponding Hamiltonian form of the Einstein equation when the boundary term does 
not vanish by using the 2+1 decomposition of the timelike S hypersurface that we have found in the f 
chapter. This of the general case will lead us to get the explicit form of the famous ADM formulas for 
conserved quantities of . 

1. 3+1 Decomposition of t lie Einstein-Hilbert Action and the corresponding 
Hamiltonian Form of the Vacuum Field Equation 

1 . 3+1 Decomposition of The Einstein-Hilbert Action 

The action for the four-dimensional vacuum field equation is of the Einstein-Hilbert action |[T9]| . ifTl 



5 = r ^R^d^x, 

J'V 



(300) 

J'V 

where the infinitesimal volume element "V is composed of the union of the neighboring hypersurfaces 
and Z,, . Symbolically, 



^ = \J^f (301) 



t=ti 
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Let us substitute the equations ( 158 1 of the 3+1 form of the spacetime Ricci scalar and ( 115 1 of the 



3+1 decomposition of g into the action ( |300| ), 
S 



R + K^ + KijK'J £rr,K DiD'N\N ^^d^x 



-i 



N[R + K'- + KijK'J] - 2£rnK - IDiD'N] ^lycTx 



(302) 



Let us convert the term JLmK into of the boundary and substitute it into the action ( 302 1 



(303) 



Then, the action (302 1 becomes 



i 

-if NVf,{Kn'')^ffd'^x. 

J'V 



N[R + K'' + KijK'J] - INV^iKn'') - INK-" - IDiD'N] ^Jyd^x 
N[R + KijK'j - K^] - IDiD'N] ^/yd'^x 



(304) 



Here we need to show that because of the boundary condition, the last integral of the action ( 304 1 
vanishes: 



f NV^,{Kni')^J^d'^x^ \ Vf,{Kn^)^P^d'^x 

J'V J'V 



(305) 



so the action ( 304 1 reduces to 



0, 



N[R + KifK'^ - K^] - 2D 



(306) 



Observe that the action (306 1 is fully composed of the intrinsic quantities of X,. This provides us to 



decompose the four-dimensional integral into of the spatial one and of the time coordinate 

-'2 



5 - ' I J [N[R + KijK'J - K^] - IDiD'N ] ^lyd'x \dt . 



(307) 



Again, the boundary term vanishes and we get the 3+1 decomposition of the Einstein-Hilbert action 
as 



(308) 
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2. The Corresponding Hamiltonian Form of the Vacuum Field Equation 



The variables of the action in the configuration we q = {■Yij,N,/3') and q = (yij,N,/3') (308 1 ifTll . 
That's, 

S = S[q,q]. 

The Lagrangian density contains the extrinsic curvature Kij [see the 3+1 decomposition of the action 



(308 1]. However, in the Hamiltonian approach, it is replaced with the configuration variables. Now, 

St 



from the 1 equation ( 184 1 of the 3+1 dimensional Einstein system, we have 



d_ 
dt 



£p\yij = -2NKi^ 



1 

- IN 



- Jij 



(309) 



We know from the ( 1 8 1 1 that X^pjij = DiPj + D j^j. Then, the equation ( 309 1 becomes 

1 



- IN 



1 

IN 



(310) 



And the Lagrangian density of the gravitational field ( 308 1 turns into 

L{q, q) = N ^(r + KjjK'^ - K^^ 



(311) 



As we see from (311), the Lagrangian density does not depend on the time derivative of N and 
so they are not dynamical variables. They are just the Lagrange multipliers. On the other hand, 
the remaining variable jij is just the dynamical variable in phase space. And the corresponding 
conjugate momentum of it is 

dL _ dKgb dL 
djij djij dKab 

= -^5ia5jbN^(y^yi^ - yj/')(6ia6jhKki + dkaSwKij) 



r'J - 



IN 

k ^ i,l ^ j k ^ i,l ^ j 



(312) 



Thus, we get 



n'j = ^\yJK-K'j\. 



(313) 
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Finally, let us find the corresponding Hamiltonian: The Legendre transformation is defined as 



(314) 



Let us substitute the explicit form of n'J (313l and y, (310l into the Hamiltonian density (314i: 



- Vy(r''^ - j( - 2NKij + Dfi + Dipj ] 
-N^Jyi^R + KijK'^ -K^^ 

- Vy{ - 2NKy^Kij + KyjDj/3i + Ky'^DiPjlNKijK'^ - K'jDj/3i 

- K'jDi/3j -NR- NKijK'j + NK^"^ 

- Vy{ - 2^K^ + KDj/3j + KDi0 + INKijK'j - K'^Dj^i 

- K'jDi/3j -NR- NKjjK'j + NK^"^ 

- Vy{ - + 2KDj/3} + NKijK'j - IK'^Dfi - A^/?| 

- Vy{ - + K^- KijK'j] + 2KDj/3j - IK^iDj^'^ 

= Vrj - n[r + k^- KijK'j] + 2[KyJi - /:^,]d/} 

- Vrj - n[R + K^- KijK'j] + 2Dj[Kyji/3' - K^i^'] 
-2p'[yjiDjK-DjKJi]] 



= Vr{ - n[r + K^- KijK'J] + 2Dj[K/3J - K^iJS'] - 2jS'[d,-^ - DjKJifj . (315) 



With the aberrations of Co = R + K^ - KjK'j and C; = DjK^i - DiK, the Hamiltonian density (315 1 
reduces to 



(316) 



<K = Vr{ - NCq + 2j3'Ci + 2Dj[K(ij - K^i/i'] 
And the related Hamiltonian is obtained by 
[ n-lcPx 



(317) 



X . 



Due to the boundary condition, the last integral of the (317 1 vanishes and we get the Hamiltonian of 
the gravitational vacuum field as 



NCQ-2p'Ci\ylyd'x, 



(318) 
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where 



Co = R + K^- KijK'i , 



(319) 



Ci = DjKji-DiK. 



(320) 



The crucial point is that the constraint equations (319\ and ^320 \ are nothing but exactly the constraint 
equations of the energy (174) and the momentum ( 176\ [for vacuum] that we deduced during the 3+1 
decomposition of the Einstein equation. Moreover, for any spacetime to be a solution of the Einstein 
equation the related Hamiltonian equation { 318\ of it must be zero. 



2. The General Gravitation Hamiltonian and The ADM Formalism 



Contrary to the previous section, we will assume that the boundary term is not zero which with the 
2+1 decomposition of the hypersurfaces that we deduced in the 1*' chapter will lead us to the well-known 
ADM formalism for the conserved quantities of a given hypersurface Ef. Here, we will assume that the 
infinitesimal four dimensional volume element d'V is the union of two spacelike hypersurfaces Xf, and Ef^ 
at the lower and upper boundaries and a timelike hypersurface S that covers the region between these two 
spacelike hypersurfaces O 

d^^^t2[ji-^t,)\Js. (321) 

Since the unit normal vector of d'V must be directed outward.However, the unit normal vector of the hy- 
persurface S,, is future-directed so it points inward. Then, with the help of the minus sign, the unit normal 
vector of 'L,^ will point outward, too. Now, in order to find the ADM formulas, we need to first deduce the 
Hamiltonian of this case: 

1 . The Gravitational Action and The Corresponding Hamiltonian when the Boundary Term is different 
than zero 

The related gravitational action is composed of the Einstein-Hilbert part, the matter (or boundary) 
part and a no dynamical part of 5 o that does not have any influence on the equation of motion H : 

SG[g]=SH[g]+SB[g]-So. (322) 

Here 



SH[g] - 




(323) 
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(324) 



(325) 



By substituting the related actions p23[ ) and ( |324[ ) into ( |322[ ), the gravitational action p22[ ) becomes 

(167r)5G = { yPgd^x + 2(£ sK ^\cPy . (326) 

Here, y" are adapted coordinates of d'V , hat are the corresponding induced 3-metric, n" are the 
corresponding unit normal vector and K is the scalar extrinsic curvature. As we mentioned before, 
because the d'V is the union of two spacehke and one timelike hypersurfaces so e = n^ria will be +1 
or -1 depending on the type of hypersurfaces. The explicit form of ( |326| ) is 

j3. 



il6n)SG= I '^R^Pgd'^x + 2 \ n"naK ^^hd^y + 2 \ n^n^KVlid' 

+ 2 f r'^ra'K^d^y 

= f ^R^Pgd'^x-2 [ K^d^y + 2 [ K^d^y 
J-v Js,2 

\ r 'K^d^y, 
JB 



y 



(327) 



+ 2 

where we have used rfria = -1 of the spacelike hypersurface and r^ra = +1 of the timelike hy- 
persurface. For convention, let us use the following form of the 3+1 decomposition of the spacetime 
Ricci scalar. 



^R = R + K"^Kab -K^- 2{n",/3n'^ - n"-//;^)^ , 



(328) 



where ";" denotes the intrinsic covariant derivative. And, we know that the 3+1 decomposition of 
the spacetime metric is 



-gd'^x = N^dtd^, 



y- 



(329) 



Then, the 3+1 decomposition of the Einstein-Hilbert part is 

J ^7? d^x = J^ ' dti^J R + K^'^Kah - \N^/hd^y 



2 f [n",pnl^ -n^nl^^p] d^x 



dt\ I R + K"^Kah-KAN^d^, 



-2^ \rf,pnl^ - . 



(330) 
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Because of d'V = U ) IJ S, we will decompose the closed integral of the equation ( 330 1 into 
the corresponding integrals of S,, , and S. For simplicity, let us work on of E,, : The spacelike 
volume element is dl,a = na^fh d^y, then, we have 



f -n"nP.a\di:, 



-2 [ K 



yfhd^ 



yfhd^ 



(331) 



where we used the fact that n"-^p is an element of the hypersurface Z^j. Similarly, by evaluating the 
corresponding integral on in the reverse direction, we will get 



; ( K^d^y. 



(332) 



Observe that the results of (|331| and (|332[) cancel out the 2"^ and 3"' terms on the right hand side of 



the gravitational action ( 332 1. Thus, the only contribution is coming from the integral over the time- 
like hypersurface S: Now, for timelike case 3-dimensional volume element is d'La - ra sf-y d^z 
and the unit spacelike normal r" of timelike hypersurface S and the unit timelike vector normal 
vector n" of the spacelike hypersurface E, are orthogonal to each other [that is, n"ra = 0], then, 

73, 



= 2 f ra-fin"i^^F^d\. 
JS 



(333) 



With the help of (333 1, the gravitational action (327 1 becomes. 



+ 2 r [K + ra■,|irfnl^]^f^d^z. 



(334) 



Notice that the the gravitational action (334i is composed of the 3+1 decomposition of the Einstein- 



Hilbert action and the integral over the timelike hypersurface S. As we mentioned before the 
Einstein-Hilbert part results in the Hamiltonian and momentum constraints that ensure whether a 
given spacetime is a solution of the The Einstein equation or not. The important point is that the 
integral over the timelike hypersurface S will lead us to the conserved quantities of the hypersur- 
faces. That's, the boundary term of S will give the ADM formalism. Therefore, we need to do the 
decompose the timelike hypersurface S by assuming that S is being foliated by the boundary of the 
spacelike hypersurface E^, St, whose topology is supposed to be S^. Now, the 2+1 decomposition 
of the 3 -metric of S is 



yd^z = Ny/(Pdtd^9. 



(335) 
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And scalar extrinsic curvature of S is 



(336) 



So with the help of ( |335[ ) and ( |335[ ), the integrand of the related integral of ( |334[ ) that is over S 
becomes 

'K + r,,pn''nl' = r„,p(g''P - r'^^) + V 



(337) 



where k is the extrinsic curvature of 5;. Thus, by substituting the result of (337 1 into the gravitational 



action (334i, we get the decomposition of the gravitational action when the boundary term is 
different than zero: 



1 

16^1 



dt{ I (R + K"''Kah-KAN^d^ 



(338) 



+ 2^ [k-kQ)N^d^eY 

Here is the extrinsic curvature of S t embedded in flat space. The is defined so that the 
gravitational action is zero for flat spacetime. 

After construction of the action, it is time to find the corresponding Hamiltonian of the system .ln 
chapter 3, we have found that the relation between the extrinsic curvature and configuration variables 

Kah ^ ^( - f'ah + Da/3b + DbPa) . (339) 

And the corresponding canonical conjugate momentum is 

d 



ab 



n 



dhab 



(340) 



Due to the fact that the boundary part is independent of hat, the equation ( 340 1 becomes 

dKmn d 



dhab 



K„ 



(341) 



where 



I6n ^Pg£G = 



R + (h"'h'"^ - h"^h"^)KahK,, 
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Then, 



2N 



■ h"''h""'Kah 



so we get 



(342) 



In order to find the corresponding Hamiltonian, let us substitute the explicit forms of tt"^ (|339[) and 



hab (342 1 in to the Hamiltonian density: 

•Kg = n"%b - ^F8£g 

-[r + K"''Kab- K^)n^ 
= yfh\lNK''''Kab - K'^Dbpa - K^'D./ib - -2NKh"^Kab 
+ Kh^^Dbpa + Kh"''Da(ib - NR - NK"^Kab + NK^ 



(343) 



^fh\NK"^Kab - NK^ -NR- 2K"''DbPa + IKh^^DbPa 
-N^[r + K^- KabK"^] -2 ^[k"^ - Kh"'']DbPa 



Co 



= -NCo - 2 ^^hlDb[fia(K"'' - Kh''^)] -pa O^^"^ - Kh"''] 



Then, by integrating the Hamiltonian density ([343 1 over the spacelike hypersurface S,, we get the 
Hamiltonian of the system as 

{\6n)HG = J \6nn-lG [k- ko)N d^O 

i^NCo-2/3aC"^^d^y (344) 
- 2 r Db[/3a{K'''' - Kh"'')] ^d^y-2(£ [k- ko)N^d^ . 
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Thus, by using the general Stokes theorem, we reach our aim of the gravitational Hamiltonian when 
the boundary term is different than zero 



-2^ [N{k - ko) + l3a{K"^ - Kh"'yb] ^/crd^O . 



(345) 



2. The ADM Formalism 
(a) The ADM Mass 

In the previous section, we have found the general gravitational Hamiltonian in the equation 
( |345 1. Due to the convention that we follow, we should do the following changes in the 



gravitational Hamiltonian (345 1 



a ^ I , b ^ J 
h ^ y , y ^ X 

k^K,ko^Ko (346) 
r ^ s , cr ^ q 

e^y. 



Then, the gravitational Hamiltonian ( 345 I turns into 



{\6n)HG = - I \NCo - 2p'Ci\^Jyd^x 



-2(£ [N{k - Ko) +l3\Kij - Ky,j)s^] ^d^y , 



(347) 



where S t is the boundary of Sf and has the topology of 5 ^ ; x is a well-defined coordinate system 
on the S/ and y is the corresponding induced 3-metric on ; a- is the scalar extrinsic curvature 
of St embedded in (S,,'y); kq is the scalar extrinsic curvature embedded in the flat spacetime 
(J.t,f)', s is the spacelike unit vector that is normal to 5r; y is a well-defined coordinate system 
on St and q is the corresponding 2-metric on Sf 

Now, suppose that a given spacetime is a solution of the Einstein equation. Then, 
the corresponding integral of constraints vanishes [ due to Co = and C, - 01 [1,1. 



H,oi. = - ^ (b A^(^ - ^o) + li\Kij - Kyij)sJ \yf^d^y. (348) 



The total mass of the which is measured by an asymptotically inertial observer (N = I and 
P -0) with a well-defined adapted coordinates of {t, x^) is given by the famous ADM energy 



78 



formula |[T1 of 



Madm ^ lim (£ (k- ko) ^(fy . 



(349) 



Madm is the conserved quantity associated the symmetry of the action under time translation 
im . And, in terms of the intrinsic connection of S, : 

Madm = Um (fi fD%- - DiifjudV ^gd''y . (350) 

Furthermore, as we mentioned before one of the conditions for a spacetime to be asymptot- 
ically flat is that there must be a coordinate system (x') in which the background metric f is 
diagonalized. Now, in this coordinate D,- = ^ and f^' = Therefore, the ADM energy 



formula (350) turns into a simpler one in this specific coordinate HJ: 

Finally, the conformal form of the ADM energy [Tl is 

Madm = -^ lim S) sHDi'V -\D^yi^^gd^y. 



(351) 



(352) 




FIG. 6. The topology of 5 .Notice that 2 is equivalent to S , and n to f in the Schwarzschild case. 

Example: The Schwarzschild spacetime in the adapted coordinates of {x") - {t, r, 6, (f>) is 

2 2 1 

g^ydxi'dx" = -(1 - -^)dt^ + (1 - dr^ + r^[de^ + sin^ddcf^] . (353) 



79 



Now, (x') = {r,6,(p) can be taken as the spatial the coordinates on Ef. Then, the induced 
3-metric is 

jij = diag[(\ - ^y\r^, ?sir?e\ . (354) 
The components of the background metric become 

9 9 9 

fij = diag{\,r ,r sin 9), (355) 

and their duals are 

= diagil, r~2, r-^sin'^e) . (356) 

It is time to start to calculate Madm- As we see in the figure ([6]), r=constant corresponds 5 f, the 
corresponding coordinates on it aiey" = {6,(p). Also, \[q(f-y = r^sin6d6d(p(dr)'. Because, 
the unit spacelike vector r is normal to the St, then, (dr)' = (1,0,0). Therefore, the related 
integral ( |350| ) of this case becomes. 



Madm - lim (f [o^yrj - Dr{f'yki)]r^sin9 dG # . (357) 
In order to evaluate the corresponding Madm that the hypersurface of the Schwarzschild holds. 



we have to first evaluate the integrands of the integral (357): Let us start to calculate 2'"^ 
integrand, 

fyu - 7rr + Xjee + 2 ^ 2^ y^0 = (l - — ) ^ + 2 . (358) 
r^sin^G \ r ' 



Since (358) is a scalar field, we have 



A(/V„) = ^(/V) = -(.-^)-^??. (359, 

Secondly, let us evaluate the 1'*' integrand: 

^^Jri = f^DkJrj = Drjrr + ^^eTrS + „ ^ -y D^r^ ■ (360) 

Now, the non-vanishing Christoffel symbols associated with D are 

fgg - -r and f^^ = -rsin^O , (361) 

Y% = fj,. - ^ and f^^ - -cosGsinG , (362) 
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With the associated covariant derivatives given by 



Dryrr = ^-2Ti,y,- = ^, (364) 
or or 



De7re = - Ty^ne - Tg^yH = -T^^yee - Tg^yrr , (365) 

D0r<l, = ^ - ^'^rTi^ri = -^^rJ^ " ^ ^^Trr , (366) 



the 1'*' integrand is 



1 1 
d ri Imx-ii-lCl n I 2m\-i 



- ^ \ (- X r'^sin^e-rsin^Oxh - —)'^\ (367) 
r^sin^O \ r ^ r ' ) 



2m/ 2m\-2 Am 1 

\ Y I f f - 2m 
2m/ 2m\-2 Am / 2m\-i 



With the help of the related results of ( |359| ) and ( |367 1, the integrand of the equation (357 1 is 
obtained from 

2m/ 2m\-2/ Am \ 
D'r,y-D,(/"r«)=-;3-(l--) + 

Am 

= — when r ^ oo . 



Let us substitute the result (368 1 into the integral equation of Madm (357 1 



Madm - hm (fi [Djyrj - Dr{f'yki)]r^smededcl> 



L r^r^sin6d0 (369) 

07T Jo Jo 



- m , 

which is exactly the mass parameter of the Schwarzschild solution. 

(b) The ADM Linear Momentum 

We have seen that one of the condition for a spacetime to be asymptotically flat is that there 
must be a Cartesian coordinate system {x') on each of Cauchy surface in which the background 
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metric is diagonal and the each diagonal element must be 1 . Therefore, the induced coordinates 
((9,)ie{ 1,2,3) for the tangent spaces Tp(L) provide three spatial direction for translation in coor- 
dinates. Now, the symmetry of the action under spatial translations provide the 2'"^ family of 
conserved quantities, !P,, of Xf. Because of the symmetry under spatial translation, we choice 
an observer for which N = and P' = I. Then, the components of related conserved quantities 
are given by the formula [ 1] of 

Pi ^^lim S (Kjk - KyjtYdiys^ ^d^y , (370) 

where the index / can take the values of (1,2,3) and it shows which component of "P,- will be 
calculated. In another words, it determines the direction in which the spatial translation will be 
done. Furthermore, the P'.s are known as the component of the ADM linear momentum of 
the hypersurface Ef of the globally-hyperbolic asymptotically flat spacetime. On the other hand, 
since the extrinsic curvature K of the embedded in the Schwarzschild spacetime with the 
standard and isotropic coordinates vanishes, then, the corresponding ADM linear momentums 
Pi of the hypersurface of the Schwarzschild spacetime vanishes. 

(c) The ADM 4-Momentum 

The ADM 4-Momentum [IJ is defined as 

K''' -(-Madm,"?), (371) 

which transform like the components of a 1-form under {x") {x ") during which the funda- 
mental properties for a spacetime to be asymptotically flat are kept valid. 

(d) The ADM Angular Momentum 

We suppose that the spacetime that we deal with receives the Killing vectors. Now, the angular 
momentum of the of a globally-hyperbolic asymptotically flat spacetime which is related to 
the rotational symmetry of the action is obtained by using the rotational Killing vectors of the 
background metric (<^;)ie{ 1,2,3) ■ In the Cartesian coordinates (x,y,z), the Killing vectors of the 
background metric about the x-axis,y-axis and z-axis are lUl 

(px = ydz - zdy , (py = zd^ - xd^ , (f>^ = xdy - ydx . (372) 

Then, the component of the angular momentum of the E, can be defined as 

Ji = ^ lim £ (Kjk - KjjuYcPiys'' ^gd^y . (373) 
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Contrary to what we expect, J'i do not transform like the 4-dimensional vectors under {x") — > 
{x ") during which the fundamental properties for a spacetime to be asymptotically flat are kept 
valid im. That's, they are coordinate-dependent. 

Because of the coordinate-dependence of J'i, the scientists have tried to put the constraints on selecting 
coordinate systems such that J'i is invariant under special subset of the related family of coordinate trans- 
formations. That's, they have being considered particular decays. For example, York |[T4]| has considered 
the following decays of the jij relative to the Cartesian coordinates for the background metric and the scalar 
extrinsic curvature 

(374) 

K = 0[r~^] . 

They are called the quasi-isotropic gauge and asymptotically maximal gauge, respectively. These asymp- 
totic gauge conditions are used to select the suitable coordinates. And it rejects some of well-known 
coordinates such as the standard Schwarzschild coordinates. Moreover, 

York proposed that the angular momentum is carried by the 0[r^^] piece of K and is invariant under the 
change within this gauge HI. 

VII. RELATION BETWEEN THE COTTON (CONFORMAL) SOLITON AND STATIC VACUUM 

SOLUTIONS 

This chapter is devoted to an application of the methods described in the previous chapters. The four- 
dimensional Einstein equation with a cosmological constant is given by 



^Rf,y - -Rg^y + Ag^,y = ^TlT . (375) 



Alternatively, we could rewrite it as 

= 8;r(r^, - i(r - ^)g,y) . (376) 



By following the same procedure as we did in chapter 4, the corresponding 3-1-1 Einstein system with A 
becomes 

(j^--Cfi)7iJ = -2NKij, (377) 



+ 4n 



A (378) 



(^-^-4;^K-25o]}> 
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R + K^ - KijK'^ ^2A + \6nE, (379) 

DjK^i - DjK - Snpi . (380) 

The static vacuum equation (y = ,p = 0) obtained from the previous system of equations is 

Rij = N~^DiDjN + Ayij, (381) 

where R = 2 A and A'^"^ Ay N = -A f20l. From now, we will use V as the intrinsic Levi-Civita connection 
and gij as the induced 3-metric (i.e. we are making a change of notation y/y — > gij and D V). The 
Cotton- York tensor lfT4l. lIBl. ifBTl is 

Cij = ej''i^,Rip-\8ip'^,R)- (382) 
Notice that the last term of the Cotton- York tensor (|382|) { M . ifBl . |[T3ll l vanishes. Let us rewrite the 



Cotton- York tensor of the static vacuum field with A in terms of the Ricci tensor: 
then 



(383) 



(384) 



NCij - ef'i[ - V,N{N-'ViVpN) + V,V,-VpA^) 

- eyP^j - [Rip - AgipWgN + VqViVpN] 

- - Agip]V,N + ejP'^V.ViVpN 

- -ejP'ilRip - AgipW.N + e/'^lLV,, V;] V^A^ + V,-V,VpA^) . 

Because the multiplication between the anti-symmettic and symmetric tensor is zero, the last term vanishes 
yielding 

NCij - -ejP'lRip - Agip]V,N + e/HV,,Vi]VpN 

- -e/n/?,> - AgipWgN - ef^R'pgiVkN (385) 

- - AgipWqN - ejPiRkpdiV'^N . 

In three dimensions the Weyl tensor is zero so the Riemann tensor can be written in terms of the Ricci 
tensor, the metric and the scalar curvature as 

Rkpqi = '2-gk[qRi]p + 2gp[iRq]k " Rgk[qgi]p ■ (386) 



By using this identity, the equation (385 1 becomes 



dj = -ep\^Rip5\ + gpiR\ - 2Agip5\]N-'\kN . (387) 



With the definition of Uu = N'^VkN, 



We can get rid off the anti- symmetric part of (388 1 : Therefore, let us first do the interchange of 



/ <^ in (388 1 



dj = -letP'iRjpUq - eif[R\ - 2K5\]Uk . 



The anti-symmetric part will drop by adding the equations ( 388 1 and ( 389 1, and we will have 



Cij = -ei'^''RjpU,-erRipU,, 



or 



Cj = -e^P^iRj^U,-erR%U,. 



Let us define Z'> = -e'PiUq and Xf = -ef^Uq, then, 



C'j = X'pRPj + R^pXf. 



This can be written as a matrix equation, let us rewrite it in the compact form 

C - XR + RX^ , 



Since C - C and R = R, the equation ( 393 1 can be written as 



C = R^X'^ + (X^)^R . 



By defining A = R and Y = X , we have 



C - A^Y + Y^A . 



Now, the matrix equation of the type 



A^X ± X'^A = B , 



has the general solution of 



1 

-I 

2 



X = ^G^BPi + G^B(1 - Pi) + (1 - P[)Y + (Pf ZP2A) , 



85 



where Z is a rank-2 antisymmetric tensor ; Pi = GA and P2 - AG such that APi - P2A - A; AGA - A 

lEa. 



By using the equation (397 1, the general solution for our equation (395 1 is 



Y = 



1, 



R C + ZR, 



or we have 



X = tCR^ - RZ. 



(398) 



(399) 



Let us examine this solution for the Cotton flow ll2T]l : The equation of the gradient Cotton soliton [21] can 
be taken as 



dj + V,-V^A^ = . 



(400) 



By choosing an ansatz, V,VyA^ = (Rij - Agij)N, the general solution for the gradient Cotton soliton ||2T1 
becomes 



2 D 
= -^-N6iP + ^{R~')iP-RrZj, 



(401) 



since X/ = -efW^ - -e/W"^ VA^, we get 



(402) 



And the corresponding constraint equation in which the rank-2 anti-symmetric tensor Z must satisfy is 

yP'^V^R/'Z^p + y^'^Ri'^V^Znp + ^-Rnnp(R'%'' - RP'Z,!) - -A . (403) 

For the case of the static vacuum solution with zero cosmological constant, we found the general solution 
as 



V"W - e'P"'Rj"Znp . 
And the corresponding constraint equation for Z is 

ejP'RfViZnp^O. 



(404) 



(405) 



It seems that the results ( |402| ), ( |403| ), ( |404| ) and ( |405| ) can be used to find which solution of the static field 
equation with A (or A = 0) is also a solution of the gradient Cotton soliton [21]. Furthermore, there is only 
one 3-dimensional Ricci soliton 11231 which is known as the Bryant soliton. However, the explicit metric 
is not known. We have not also been able to solve the constraint equations and have not found explicit 
metrics. But the formulation outlined above can be used to explore the gradient Cotton solitons f2T\ and 
the solutions of Topologically Massive Gravity (TMG) Il24ll . ||25l . |[26l as well as the gradient Ricci solitons 
1123. 
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VIII. CONCLUSION 



In this work, we have first learned how to foliate a globally hyperbolic four-dimensional spacetime 
by a continuous set of Cauchy surfaces (Lt)te'R which, with particular numbers of the projection onto the 
hypersurface and along the unit normal vector, provide us the fundamental relations of the 3+1 formalism 
|[T2l . Furthermore, with the help of these basic relations as well as the 3+1 decomposition of the stress- 
energy tensor, we have learned that the full projection of the Einstein equation onto X gives the dynamical 
part and the other two projection give the constraint equations which are used to check whether a given 
spacetime is a solution of the Einstein equation or not. Moreover, with the help of the coordinate adapted to 
the flow and the shift vector yS, one can convert the 3+1 Einstein system to a set of PDEs lU. 

Secondly, we have analyzed the flow of the hypersurfaces as if there is a conformal relation between a 
well-defined conformal background metric y and the set of the induced 3-metrics associated with the hyper- 
surfaces. Moreover, by constructing the fundamental conformal transformations of the intrinsic quantities 
of the hypersurface, the 3+1 conformal Einstein system is constructed. And we have emphasized that the 
trace and traceless parts of the 3+1 dynamical Einstein equation transform separately under the conformal 
transformation. Finally, we have reconstructed the 3+1 Einstein system for the foliation which is maxi- 
mally sliced (K = 0). We have seen that the conformal background metric is nothing but a conformally flat 
background metric and the Cotton- York tensor |[T3l . lfT4l vanishes in this case. We have also seen that this 
particular case leads to the 3+1 IWM system iflTl . ifTSi which is the conformal approximation to the general 
relativity. 

Thirdly, we have reconstructed the Hamiltonian form ||9]|, lITOll of the general relativity which provides 
us the ADM formalism for the conserved quantity of hypersurfaces of the globally-hyperbolic spacetimes 
which asymptotically approach to the well-defined spacetimes such as the Minkowski spacetime lfT6l . Fur- 
thermore, we have seen that the quasi-isotropic gauge and the asymptotically maximal gauge force us to 
shrink the cluster of the coordinates in which the ADM angular momentum J'l becomes invariant. 

Finally, we have proposed a method in chapter 7 which we think will give the relation between the 
solutions of the gradient Cotton soliton ll2T1l and of the static vacuum field equations. Furthermore, we 
think that this method can be used to find the relation between the solutions of the gradient Ricci ||23]| and 
the Cotton ini solitons and the solutions of the Topologically Massive Gravity (TMG) 1211, ||25]|, HH. 
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